NAVAL K 
MONTEREY. 



-m-'/ j tgRARY 

i- JE SCHOOL 
OrtMlA 93943 



I 

NAW 

MOf 



NFS-69-84-008 



NAVAL 



POSTGRADUATE SCHOOL 

Monterey, California 




FILMWISE CONDENSATION OF STEAM ON 
LOW INTEGRAL- FINNED TUBES 

by 

Ioannis V. Georgiadis 
September 1984 



Thesis Advisor: P. J. Marto 



Approved for public release; distribution unlimited 
Prepared for: 

National Science Foundation 
Division of Engineering 
Washington, DC 20550 



T22 




NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



Commodore Robert H. Shumaker David A. Schrady 

Superintendent Provost 



This thesis prepared in conjunction with research 
supported in part by National Science Foundation, Division 
of Engineering, Washington, DC under MEA82-03567. 

Reproduction of all or part of this report is authorized. 



Released as a 
Technical Report by: 



UNCLASSIFIED 



SECURITY CLASSIFICATION of this PAGE (When Data Entered) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


t. REPORT NUM8ER 

NPS-69-84-00 8 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


4. TITLE (and Subtitle) 

Filmwise Condensation of Steam on Low 
Integral-Finned Tubes 


5. TYPE OF REPORT 4 PERIOD COVERED 

Engineer's Thesis 
September 1984 


6. PERFORMING ORG. REPORT NUMBER 


7. AU THORf a) 

Ioannis V. Georgiadis 

s 


8. CONTRACT OR GRANT NUMBERfsJ 

Agreement No. MEA82-03567 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93943 


10. PROGRAM ELEMENT, PROJECT, TASK 
AREA 4 WORK UNIT NUMBERS 


11. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93943 


12. REPORT DATE 

September 1984 


13. NUMBER OF PAGES 
166 


U. MONITORING AGENCY NAME & ADDRESS^/ different from Controlling Office) 

National Science Foundation 
Washington, dc 20550 


is. SECURITY CLASS, (ot thte report) 

Unclassified 


IS*. DECLASSI F| CATION / DOWN GRADING 
SCHEDULE 



16. DISTRIBUTION STATEMENT (of this Report) 



Approved for public release; distribution unlimited 



17. DISTRIBUTION STATEMENT (of the abstract entered In Block 20, If different from Report) 



18. SUPPLEMENTARY NOTES 



19. KEY WORDS (Continue on reverse aide if neceeaary and identify by block number) 

Steam, Condensation, Filmwise, Tube, Low-Integral Fin, 
Heat-Transfer Coefficient, Enhancement 



20. ABSTRACT (Continue on reverse aide If neceaeary and identify by block number) 

Filmwise condensation heat-transfer measurements of steam 
were made on horizontal tubes under vacuum and near-atmospheric 
pressures. Data were taken for a smooth tube and for 21 tubes 
which contained rectangularly-shaped, low integral fins. The 
fin geometry was systematically varied in order to investigate 
the dependence of the steam-side heat-transfer coefficient on 
fin spacing, thickness and height. 



UNCLASSIFIED 



DD 1 JAN 73 1473 EDITION OF 1 NOV 65 IS OBSOLETE 

S N 0102- LF- 014- 6601 1 



SECURITY CLASSIFICATION OF THIS PACE (When Data Bntarad) 



SECURITY CLASSIFICATION OF THIS PAGE (Whft Data Enffd) 



#20 - ABSTRACT - (CONTINUED) 



The condensation process was found to be most sensitive 
to fin spacing and to be weakly dependent on fin thickness. 
When the fin height was increased from 1.0 to 2.0 mm, the 
increase in the heat-transf er coefficient was about 20 
percent, while the increase in condensing area was 51 
percent . 

The best performing finned tube had a fin spacing of 1.5 
mm, a fin thickness of 1.0 mm and a fin height of 2.0 mm. 

It showed an enhancement over the smooth tube of between 4 
to 5 under vacuum conditions and around 6 at atmospheric 
conditions. Thus, the use of finned tubes may result in a 
significant reduction in the size and weight of present-day 
steam condensers. 



S' N 0102- LF- 014- 6601 



UNCLASSIFIED 



2 



SECURITY CLASSIFICATION OF THIS PAGEf»?i«n Data Entarad) 



Approved for public release; distribution unlimited. 



Filmwise Condensation of Steam on 
Low Integral-Finned Tubes 

by 

Ioannis V. georgiadis 
Lieutenant, Hellenic Navy 
B.S., Hellenic Naval Academy, 1975 



Submitted in partial fulfillment of the 
requirements for the degrees of 



MASTER OF SCIENCE IN MECHANICAL ENGINEERING 

and 

MECHANICAL ENGINEER 
from the 

NAVAL POSTGRADUATE SCHOOL 
September 1984 



»„S.gS“»T 

N MOttTEREY, CALIFORNIA 93943 



(s 7 £ 5 -L 

v • ' 



Filmwise condensation heat- transfer measurements of 
steam were made on horizontal tubes under vacuum and near- 
atmospheric pressures. Data were taken for a smooth tube 
and for 21 tubes which contained rectangularly-shaped, low 
integral fins. The fin geometry was systematically varied 
in order to investigate the dependence of the steam-side 
heat- transfer coefficient on fin spacing, thickness and 
height . 

The condensation process was found to be most sensitive 
to fin spacing and to be weakly dependent on fin thickness. 
When the fin height was increased from 1.0 to 2.0 mm, the 
increase in the heat-transfer coefficient was about 20 
percent, while the increase in condensing area was 51 
percent . 

The best performing finned tube had a fin spacing of 1.5 
mm, a fin thickness of 1.0 mm and a fin height of 2.0 mm. 

It showed an enhancement over the smooth tube of between 4 
to 5 under vacuum conditions and around 6 at atmospheric 
conditions. Thus, the use of finned tubes may result in a 
significant reduction in the size and weight of present-day 
steam condensers. 



4 



TABLE OF CONTENTS 



MOW 



1 iOiL i , C/-»Ui Ou../a 






n 

u 



I. INTRODUCTION 17 

A. FILMWISE CONDENSATION ON FINNED TUBES 18 

B. SURVEY OF ANALYTICAL MODELS AND 

EXPERIMENTAL WORK 22 

C. OBJECTIVES 28 

II. DESCRIPTION OF APPARATUS 30 

A. SYSTEM OVERVIEW 30 

B. SYSTEM INSTRUMENTATION 33 

C. DATA ACQUISITION SYSTEM 34 

D. SYSTEM MODIFICATIONS 35 

E. VACUUM INTEGRITY 38 

F. TUBES TESTED 39 

G. INSERT 42 

III. SYSTEM OPERATION AND PROBLEMS 45 

A. THE DROPWISE CONDENSATION PROBLEM 47 

B. VAPOR VELOCITY LIMITATIONS 49 

IV. DATA REDUCTION 51 

A. PROGRAM SIEDER 53 

B. MODIFIED WILSON PLOT PROGRAM 55 

C. MAIN DATA REDUCTION PROGRAM ( DRP4 ) 56 

V. RESULTS AND DISCUSSION 57 

A. INSIDE HEAT-TRANSFER COEFFICIENT 57 

1. Instrumented Tube Results 57 

2. Uninstrumented Tube Results 64 

3. Summary of the Inside Heat-Transfer 

Coefficient 66 



5 



B. PERFORMANCE OF FINNED TUBES 69 

1. Condensation on Fins in General 69 

2. Significance of Fin Spacing 75 

3. Significance of Fin Thickness 99 

4. Significance of Fin Height — 107 

5. Significance of Tube Diameter 116 

VI. CONCLUSIONS AND RECOMMENDATIONS 119 

A. CONCLUSIONS 119 

B. RECOMMENDATIONS 120 

APPENDIX A: PREPARATION OF CONDENSER TUBES 122 

APPENDIX B: SYSTEM START-UP AND SHUT-DOWN PROCEDURES - 123 

APPENDIX C: RAW DATA OF UNINSTRUMENTED AND 

INSTRUMENTED TUBE 125 

APPENDIX D: RAW DATA OF TUBES TESTED 127 

APPENDIX E: ERROR ANALYSIS . 150 

LIST OF REFERENCES 163 

INITIAL DISTRIBUTION LIST 165 



6 



LIST OF TABLES 



I. Geometry of Finned Tubes 42 

II. Instrumented Tube Results for the Vacuum Run 61 

III. Instrumented Tube Results for the 

Atmospheric Run 7 61 

IV. Constants of Equation (5.6) for All Tested 

Tubes 7 2 

V. Raw Data of Uninstrumented Smooth Tube 125 

VI. Raw Data of Instrumented Tube 126 

VII. Raw Data of Tube with Fins of s = 0.5, t = 1.0 

and e = 1.0 mm 12 8 

VIII. Raw Data of Tube with Fins of s = 1.0, t = 1.0 

and e = 1.0 mm 129 

IX. Raw Data of Tube with Fins of s = 1.5, t = 1.0 

and e = 1.0 mm 13 0 

X. Raw Data of Tube with Fins of s = 2.0, t = 1.0 

and e = 1.0 mm 131 

XI. Raw Data of Tube with Fins of s = 4.0, t = 1.0 

and e = 1.0 mm 13 2 

XII. Raw Data of Tube with Fins of s = 0.5, t = 0.75 

and e = 1.0 mm 133 

XIII. Raw Data of Tube with Fins of s = 1.0, t = 0.75 

and e = 1.0 mm 134 

XIV. Raw Data of Tube with Fins of s = 1.5, t = 0.75 

and e = 1.0 mm 135 

XV. Raw Data of Tube with Fins of s = 2.0, t = 0.75 

and e = 1.0 mm 136 

XVI. Raw Data of Tube with Fins of s = 4.0, t = 0.75 

and e = 1.0 mm 137 

XVII. Raw Data of Tube with Fins of s = 0.5, 7 = 0.5 

and e = 1.0 mm 138 



7 



XVIII . 

XIX. 

XX. 

XXI. 

XXII. 
XXIII. 

XXIV. 

XXV. 

XXVI. 
XXVII. 
XXVIII. 



Raw Data of Tube with Fins of s = 1.0, t = 0.5 
and e = 1.0 mm 139 

Raw Data of Tube with Fins of s = 1.5, t = 0.5 
and e = 1.0 mm 140 

Raw Data of Tube with Fins of s = 2.0, t = 0.5 
and e = 1.0 mm 141 

Raw Data of Tube with Fins of s = 4.0, t = 0.5 
and e = 1.0 mm 142 

Raw Data of Tube with Fins of s = 1.0, t = 1.0 
and e = 2.0 mm 14 3 

Raw Data of Tube with Fins of s = 1.5, t = 1.0 
and e = 2.0 mm 144 

Raw Data of Tube with Fins of s = 2.0, t = 1.0 
and e = 2 . 0 mm 145 

Raw Data of Tube with Fins of s = 4.0, t= 1.0 
and e = 2 . 0 mm 14 6 

Raw Data of Tube with Fins of s = 1.0, t = 1.5 
and e = 1 . 0 mm 14 7 

Raw Data of Tube with Fins of s = 1.25, 
t = 1.2 5 and e = 1.0 mm 148 

Raw Data of Smooth Tube 149 



8 



LIST OF FIGURES 



1.1 Condensate Profile at Top Part of Tube 19 

1.2 Typical Schematic of Retention Angle 24 

2.1 Schematic of Test Apparatus 31 

2.2 Schematic Details of Test Section (insert 

not Shown) 3 2 

2.3 Schematic Details of Purge System and Water 

Sump 3 7 

2.4 Photograph of the Instrumented Tube 40 

2.5 Photograph of a Set of Tested Tubes with 

Fins of t = 0.75 and e = 1.0 mm 41 

2.6 Photograph of the Insert Together with a 

Finned Tube Tested 44 

4.1 Mixing Chamber Calibration Curve 52 

5.1 Temperature Distribution for Instrumented 

Tube (Vacuum Run) 59 

5.2 Temperature Distribution for Instrumented 

Tube (Atmospheric Run) 60 

5.3 Correlation of the Water-Side Coefficient 

Obtained with Instrumented Smooth Tube 63 

5.4 Comparison of Smooth-Tube Data with Fujii 

Correlation 67 

5.5 Variation of Heat Flux with Steam-Side 

Temperature Drop 70 

5.6 Variation of Condensing Coefficient with 

Heat Flux 73 

5.7 Cross Plot Showing Best Fin Spacing Among 

Tubes with t = 1.0 and e = 1.0 mm 77 

5.8 Miniscus Profile for Various Fin Sapcings 

(a) s = 1.0 mm, (b) s = 1.5 mm, (c) s = 2.0 mm -- 79 



9 



5.9 Normalized Enhancement Ratio Versus Fin 
Spacing for Tubes with Fins of t = 1.0 mm 

and e = 1 . 0 mm 80 

5.10 Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 1.0 

and e = 1.0 mm (Vac. Runs) * 83 

5.11 Variation of Heat Flux with Steam-Side 
Temperature Drop for Tubes with t = 1.0 and 

e = 1.0 mm (Vac. Runs) 84 

5.12 Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 1.0 

and e = 1.0 mm (Atm. Runs) 85 

5.13 Variation of Heat Flux with Steam-Side 
Temperature Drop for the Set of Tubes with 

t = 1.0 and e = 1.0 mm (Atm. Runs) 86 

5.14 Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 0.75 

and e = 1 . 0 mm (Vac. Runs) 87 

5.15 Variation of Heat Flux with Steam-Side 
Temperature Drop for the Set of Tubes with 

t = 0.7 5 and e = 1.0 mm (Vac. Runs) 8 8 

O 

5.16 Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 0.75 

and e = 1.0 mm (Atm. Runs) 89 

5.17 Variation of Heat Flux with Steam-Side 
Temperature Drop for the Set of Tubes with 

t = 0.75 and e = 1.0 mm (Atm. Runs) 9 0 

5.18 Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 0.5 

and e = 1.0 mm (Vac. Runs) 91 

5.19 Variation of Heat Flux with Steam-Side 
Temperature Drop for the Set of Tubes with 

t = 0.5 and e = 1.0 mm (Vac. Runs) 9 2 

5.20 Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 0.5 

and e = 1.0 mm (Atm. Runs) 9 3 

5.21 Variation of Heat Flux with Steam-Side 
Temperature Drop for the Set of Tubes with 

t = 0.5 and e = 1.0 mm (Atm. Runs) 9 4 



10 



5.22 Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 1.0 

and e = 2.0 mm (Vac. Runs) 95 

5.23 Variation of Heat Flux with Steam-Side 
Temperature Drop for the Set of Tubes with 

t = 1.0 and e = 2.0 mm (Vac. Runs) 9 6 

5.24 Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 1.0 

and e = 2.0 mm (Atm. Runs) 97 

5.25 Variation of Condensing Coefficient with 
Heat Flux for the Set of Tubes with t = 1.0 

and e = 2 . 0 mm (Atm. Runs) 98 

5.26 Cross Plots Showing Best Fin Spacing Among 
Sets of Tubes with Fins of t = 0.5, 0.75 and 

1 . 0 mm 10 0 

5.27 Influence of Fin Thickness on Enhancement 
Ratio for Tubes with s = 1.0 mm, for both 

Pressure Conditions 102 

5.28 Variation of Condensing Coefficient with 
Heat Flux for Tubes with Fin Geometry Shown 

in the Figure (Vac. Runs) 10 8 

5.29 Variation of Condensing Coefficient with 
Heat Flux for Tubes with Fin Geometry Shown 

in the Figure (Atm. Runs) — 109 

5.30 Cross Plot Showing Best Fin Spacing for 
Tubes with e = 2.0 mm and t = 1.0 mm, for 

both Pressure Conditions 110 

5.31 Comparison of Enhancement Ratio Obtained on 
Sets of Tubes with e = 1.0 and e = 2.0 mm 

(Vac. Runs) HI 

5.32 Condensate Profiles for Tubes with Fins of 

Different Height 113 

5.33 Normalized Enhancement Ratio of Sets of 

Tubes with Fins of 1.0 and 2.0 mm Height 115 

5.34 Cross Plots Showing Best Fin Spacing for 

Tubes of Different Diameter 117 

E.l Error Bands for a Tube Tested Under Vacuum 

Conditions 153 

E.2 Error Bands for a Tube Tested Under 

Atmospheric Conditions 154 



11 



LIST OF SYMBOLS 




e 

F 



g 



h 



BK 



h 



c 



h 



eff 




hi 



Effective Area of Finned Tube (See Eqn. 1.10) 
Outside Area of Finned Tube 
Outside Area of Smooth Tube 

Constant Used with Sieder-Tate- type Equation 

Constant Used with Sieder-Tate- type Equation 

Correction Factor for Cooling Water Viscosity 
Changes 

Constant Pressure Specific Heat of Cooling Water 

Equivalent Diameter Used with Beatty and Katz 
Equation 

Diameter to Outside of Fins 

Inside Diameter of Condenser Tube 

Outside Diameter of Smooth Tube, Equal to Root 
Diameter of Finned Tube 

Fin Height 

gD y.eh.c /V^k^AT 
^ of fg vf 

Local Gravitational Acceleration 

Average Heat- transfer Coefficient Given by 
Beatty and Katz Model 

Heat-transfer Coefficient Based on Smooth Tube 
Diameter 

Heat- transfer Coefficient on the Flooded part 
of the Tube 

Specific Enthalpy of Vaporization 

Water-side Heat- transfer Coefficient 

Effective Heat- transfer Coefficient on the 
Flooded Part of the Tube 



12 



h ,ho 
o 


Steam-side Heat-transfer Coefficient 


^ow 


Average Heat-transfer Coefficient Given by 
Owen et al. Model 


h FW 


Average Heat- transfer Coefficient Given by Rudy 
and Webb Model 


h 

u 


Heat-transf er Coefficient on the Upper Tube 
Portion 


k 

c 


Thermal Conductivity of Cooling Water 


k eff 


Effective Thermal Conductivity 


k f 


Thermal Conductivity of Condensate 


k fin 


Thermal Conductivity of Fin Material 


m act 


Actual Cooling Water Mass Flow Rate 


m calc 


Calculated Cooling Water Mass Flow Rate 


L 


Condensing Tube Length 


LMTD 


Logarithmic-mean- temperature Difference 


N 


Fin Density 


Nu 


Water-side Nusselt Number, h . D . /k 

11C 


Nu 

g 


Steam-side Nusselt Number, h D /k- 

oof 


ap ab' ap cb 


Pressure Difference Between Points C and B 
or Points A and B (Fig. 1.1) 


P A ,P B ,P C 


Pressure in Condensate (Fig. 1.1) 


Pr 


Prandtl Number of Cooling Water, C u /k 

^ pc H c 7 c 


P 

V 


Pressure of Steam 


q 


Heat Flux Based on Outside Area of Smooth Tube 


Re 


Coolant-side Reynolds Number, P c v c D -j_/^ c 


Re 


Steam-side Reynolds Number, p^V D o /y£ 


s 


Fin Spacing 


t 


Fin Thickness 



13 



Tin 

T or Ts 
s 

Tout 



T 



wo 



AT 

AT 




U 



V 

V 
X 



O 



v 



Vw 



x 



y 



Cooling Water Inlet Temperature 

Temperature of Steam 

Cooling water Outlet Temperature 

Average Outer all Temperature, Temperature at 
Root of Fins 

Local Temperature Drop Across Condensate Film 

Average Temperature Drop Across Condensate Film 

(T - T ) 
s wo 

Cooling Water Temperature Drop 

Radius of Curvature of the Miniscus Profile in 
the Interfin Space at the Top Part of the Tube 
(Fig. 1.1) 

Overall Heat-transfer Coefficient 
Cooling Water Velocity 
Steam Velocity 

Average Fin Height (Defined by Equation (1.11)) 
Coordinate Axis (Fig. 1.1) 

Coordinate Axis (Fig. 1.1) 



Greek symbols 

0 Angle Measured from Top of Tube 

y Viscosity of Cooling Water at Bulk Temperature 

c 

y^ Viscosity of Condensate at Film Temperature 

y Viscosity of Cooling Water at Inner Wall 

Temperature 

p Density of Cooling Water 

Pf Density of Condensate 

P v Density of Steam 

o Surface Tension of Condensate 



14 



Fin Efficiency 



ip Condensate Retention Angle Measured from the 

Bottom of the Tube to the Position at Which 
the condensate Just Fills the Interfin Space 



15 



ACKNOWLEDGEMENTS 



This thesis was completed because of the careful and 
good-willing help of the following persons : 



1 . 


Dr. 


P. 


J. 


Marto , my advisor. 


2. 


Dr. 


J. 


W. 


Rose, co-advisor in the project. 


3. 


Dr. 


A. 


S . 


Wanniarachchi , my second advisor 



3. Mr. Tom Christian and the workers in the machine 
shop of the Mechanical Engineering Department, and 

4. Artemis, my wife. 

I would like to express to them my deep feelings by just 
saying "Thank you very much." 



16 



I . 



INTRODUCTION 



Modern marine propulsion plants must be designed to 
provide maximum power with small size, low weight and 
reasonable cost. A heavy power plant may have an impact on 
ship dimensions and displacement since the power required to 
achieve a given speed increases with the displacement. A 
sizeable propulsion plant will occupy space that is valua- 
ble, while the machinery room will have poor accessibility. 

A significant decrease in size and weight of a steam 
power plant can be achieved by reducing the size of the 
condenser. Present condensers are large because of the use 
of smooth tubes and generous design margins. The heat- 
transfer performance of any condenser tube is limited by the 
thermal resistances on the water-side, in the wall and on 
the steam-side of the tube. Reduction of any of these 
resistances will contribute to a lighter and smaller 
condenser . 

In recent years, various efforts have been made in order 
to increase the outside (i.e., steam-side) heat- transfer 
coefficient, and these methods have been thoroughly reviewed 
by Marto [Ref. 1] . Outside enhancement techniques include 
the use of low integral fins, roped tubes, fluted tubes, 
drainage strips attached on the tubes and applied coatings 
to promote dropwise condensation. 
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A. FILMWISE CONDENSATION ON FINNED TUBES 

Filmwise condensation heat transfer on a smooth tube is 
limited by the thickness of the condensate layer. This 
layer forms a relatively high conduction resistance between 
the vapor and the tube surface. To enhance filmwise conden- 
sation heat transfer, it is then important to thin this 
condensate film. One common way to decrease this condensate 
film is to use radial fins attached to the external surface 
of the tube in order to promote surface-tension effects. 

The theoretical treatment of this problem is very diffi- 
cult owing to the large number of controlling parameters 
such as the three-dimensional condensate flow pattern, vapor 
shear, surface-tension forces, wall conduction effects and 
condensate retention. Studies following Gregorig [Ref. 2] 
have demonstrated the important role of surface tension in 
modifying the condensate film over the fin surface. The 
presence of fins causes three main effects: 

1. the condensing area increases, 

2. the surface-tension forces thin the condensate film 
over the tip of the fins causing improved heat 
transfer, and 

3. the fins may cause condensate retention especially 
on the lower part of the tube. 

The first two effects are beneficial, while the third 
undermines the heat-transfer performance. The second effect 
is most important in the upper portion of the tube, where 
there is little condensate retention. For steam condensa- 
tion data on finned tubes, Yau et al. [Ref. 3] measured 
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enhancement ratios (i.e., the ratio of vapor-side coeffi- 
cient of a finned tube to the value of smooth tube) in 
excess of the increase in condensing area. This additional 
enhancement is most probably due to the thinning effect of 
the surface-tension forces acting on the condensate film. 
This phenomenon can be explained by referring to Figure 1.1. 



Condensate 




Figure 1.1. Condensate Profile at Top Part of Tube 



Since the condensate film has a concave shape at point 
A, the pressure within the condensate at this point is 
smaller than the local vapor pressure, and is given by: 



A similar expression may 



P v - r 



a 



( 1 . 1 ) 



be written at point C 



P 



c 




( 1 . 2 ) 
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Note that the local radii of curvature at points A and C are 
very large; therefore, the liquid pressure at these points 
nearly equals the vapor pressure. On the other hand, at 
point B the condensate profile has a more concave shape, 
thus resulting in a lower pressure than at either points A 
or C : 



P 



B 




( 1 . 3 ) 



As a result, there exists a pressure difference from point A 
to point B as given by: 



AP 



AB 




( 1 . 4 ) 



and in a similar way from point C to point B as given by: 



AP 



CB 




( 1 . 5 ) 



In general, these pressure differences are positive, causing 
the condensate to be pulled toward point B (which acts as a 
condensate run-off channel) , thus thinning the condensate 
film on the fins as well as in the interfin area. This 
thinning effect, therefore, increases the heat-transfer 
performance in excess of the surface area gain. 
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On the other hand, as mentioned above, the third effect 
of the fins which is to retain the condensate, especially on 
the lower part of the tube, creates a relatively thick 
condensate film, causing lower heat transfer. When the fin 
spacing is reduced, a high surface-tension condensate will 
bridge the fins; therefore, the heat- transfer performance 
may decrease despite that the condensing area has increased. 
Thus, owing to the competing nature of these effects, the 
selection of a best fin geometry needs considerable attention. 

The effect of fin spacing and fluid properties (such as 
surface tension) on condensate retention has been Well 
studied by various investigators [Refs. 4,5,6]. A detailed 
discussion of these studies is presented in the next section. 

One more important effect that is common for both smooth 
and finned tubes is the oscillation of the condensate layer 
resulting from the intermittent departure of condensate 
droplets from the lower part of the tube surface. When a 
droplet departs, the nearby film stretches (becomes thinner) 
and higher heat transfer takes place locally. After the 
droplet has departed, the film grows again and becomes 
thicker. The frequency and amplitude of this oscillation 
are functions of the heat flux, fin geometry and the proper- 
ties of the fluid. To the best of this author's knowledge, 
this transient phenomenon has not yet been investigated. 
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B. SURVEY OF ANALYTICAL MODELS AND EXPERIMENTAL WORK 

In 1948, Beatty and Katz [Ref. 7] developed a simple 
analytical model to predict the condensing heat-transfer 
coefficient on a horizontal finned tube. Their model was 
supported by various other researchers, as stated in 

[Ref. 8] . The final form of the Beatty and Katz prediction 

% 

has the following equations: 



h 



Bk 



= 0.689 



j^p 



I 



f Kf( Pf -Pylghfg 

y f at 




( 1 . 6 ) 



where 



1/4 



e 



A ef f D 1/M 
o 



4>A 4 



+ 1.3 



A ef f X 1/4 



(1.7) 



A = 7T D L 
r o 



( 1 . 8 ) 



A 



f 




d 2 )nl 

o 



(1.9) 



A r r = A + <j)A_, and 
ef f r f 



X 




4 D 



f 



( 1 . 10 ) 



( 1 . 11 ) 



The major fault with their model was the neglect of any 
surface-tension effects. They therefore assumed that the 
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Nusselt equation can be applied to all the external surfaces 
of the tube (i.e., tip, sides of fins and interfin space). 
According to [Ref. 8 and 9], this model overpredicts the 
condensing coefficient especially for those fluids that have 
high surface tensions (for example, water). 

Rudy et al. [Ref. 5] attempted to include the surface- 
tension effect by modifying the original Beatty and Katz 
equation. To obtain this modification, they first investi- 
gated the retention angle problem. They performed static 
(i.e., with no condensation) and dynamic (i.e., with conden- 
sation) measurements of the retention angle on tubes with 
fin densities from 0.75 to 1.38 fins/mm. They found that 
the retention angle (see Figure 1.2) can be very well pre- 
dicted (within ±10 percent) by the following equation: 
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The above equation was tested against many fluids such as 
refrigerant-11, pentane and water. It is worth mentioning 
that Honda et al. [Ref. 6] and Owen et al. [Ref. 8], 
obtained essentially the same equation by using different 
models. Next, Rudy et al. modelled the condensation process 
assuming no heat transfer in the flooded portion of the 
tube. They developed the following prediction: 
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where h D „ is given by equation (1.6). Unfortunately, this 
model tends to underpredict most of the available experimental 
data by as much as about 60 percent [Ref. 9]. 

Owen et al. [Ref. 8] modified the Rudy et al. model by 
including heat transfer from the flooded part of the tube. 

For the fluids that they tested, their model predicted the 

% 

condensing coefficient within ±30 percent. The prediction 
that they developed is based on the following equations : 
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The above analysis uses Beatty and Katz prediction (equation 
(1.6)) for the unflooded part of the tube. As mentioned 
earlier, the surface-tension effect is not included in equa- 
tion (1.6). However, as discussed in the previous section 
(Section A) , the surface-tension forces play an important 
role in filmwise condensation by thinning the condensate 
film. Therefore, the Owen model needs further improvements 
to account for the surface-tension effect at the unflooded 
part of the tube. 

Recently, Honda et al. [Ref. 9], developed a model to 
predict -the condensing coefficient by accounting for surface- 
tension effects and for the non-isothermal behavior of the 
fins. Their model can predict most of the available experi- 
mental data including 11 fluids and 22 tubes within ±20 
percent. Honda et al . provide a comparison of these four 
models with experimental data. It is worth mentioning that 
the worst disagreement (within ±30 percent) between the data 
and their model occurs for the steam data of Rose [Ref. 3]. 
This is probably owing to the high surface tension of water, 
as the surface-tension effects are still inadequately 
accounted for even in this model. Further, it is shown that 
all the models except that of Honda overpredict these steam 
data. In comparison with the other models, the Honda model 
shows considerable promise. However, it must be noted that 
this model is rather laborious as it requires considerable 
computer time for iterative finite difference techniques. 



26 



The Honda model may require further improvements to predict 
steam data more accurately. Also, most importantly, a more 
practical model must be developed so that the designers may 
use it with confidence. Such a model may consist of a' 
combination of a few correlations. Thus, it is extremely 
important to generate a sufficient amount of experimental 
data covering all important geometric variables (such as fin 
spacing, fin thickness, fin height, fin shape, tube diameter, 
etc.) for condensation of steam on finned tubes. 

In 1961, Staub [Ref. 10] conducted steam condensation 
experiments with horizontal fluted tubes under vacuum and 
atmospheric conditions. He used a bell-shaped (rather than 
rectangular) fin profile on a root diameter of 14.9 mm. These 
fins had a pitch of 1.0 mm, a thickness of 0.4 mm at the 
root of the fin and a height of 0.5 mm. He reported a vapor- 
side heat- transfer-coefficient enhancement (i.e., the ratio of 
vapor-side coefficient of finned tube to the value of smooth 
tube) of around 3 for both vacuum and atmospheric pressure 
conditions . 

At Queen Mary College of London, Yau et al. [Ref. 3] obtained 
steam condensation data at atmospheric pressure on various 
rectangularly-shaped finned tubes having a root diameter of 
12.7 mm. They found an optimum spacing of 1.5 mm for a fin 
thickness of 0.5 mm and a fin height of 1 mm. Their optimum 
tube resulted in a vapor-side heat- transfer coefficient 
enhancement of 3 to 4 . 

On-going research is underway at the Naval Postgraduate 
School to study the feasibility of using low integral-finned 



27 



tubes with steam. A test apparatus has been constructed by 
Krohn [Ref. 11] to systematically study steam condensation 
on a single horizontal tube. The test apparatus was instru- 
mented and tested by Graber [Ref. 12], and further system 
improvements for leak tightness were made by Poole [Ref. 13]. 
Poole conducted experiments both under vacuum and at pressures 
near-atmospheric. He used machined, rectangular fins with 
a fin thickness and a fin height of 1 mm on a root diameter 
of 19 mm. He tested a total of six finned tubes with fin 
spacings of 0.5, 1.0, 1.5, 2.0 , 4.0 and 9.0 mm as well as a 
smooth tube. His results showed that the best enhancement 
of 3 to 4 occurred for the tube with a fin spacing of 1.5 mm. 
However, Poole had problems with partial dropwise condensa- 
tion since the total filmwise mode lasted only for about 20 
minutes during his experimental runs. Thus, it became 
necessary to repeat his test runs in order to verify his 
preliminary results. 

C. OBJECTIVES 

The main objectives of this research effort were, there- 
fore, to: 

1) establish a preparation procedure that will ensure 
filmwise condensation on copper finned tubes. 

2) take data on an instrumented smooth tube to determine 
the water-side heat-transfer coefficient with internal 
enhancement . 
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3) take data on a smooth tube to check the data-reduction 
program and to establish a basis for comparison with the 
finned tubes. 

4) take data on tubes with rectangularly shaped fins of 
various fin spacings, fin thicknesses and fin heights to 
obtain a possible optimum fin geometry. 
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II. DESCRIPTION OF APPARATUS 



A. SYSTEM OVERVIEW 

The apparatus used for this research was essentially the 
same used in [Ref. 13], with several small modifications. A 
schematic of the system is shown in Figure 2.1. Steam was 
generated in a 304.8 mm (12 in) diameter Pyrex glass boiler 
using ten 4000 Watt, 480 Volt Watlow immersion heaters. 
Passing through a 304.8 mm (12 in) to 152.4 mm (6 in) 
reducing section, the steam travelled upward through a Pyrex 
glass section 2.44 m (8 ft) in length, around a 180-degree 
bend and back down a straightening section 1.52 m (5 ft) in 
length before entering the stainless-steel test section. 

The condenser tube to be tested was mounted horizontally in 
the test section behind a viewport to permit visual observa- 
tion of the condensing process. Figure 2.2 shows the 
details of the test tube mounted in the test section. 

Steam that did not condense on the test tube passed into 
a stainless steel auxiliary condenser, and all condensate 
was returned via gravity to the boiler. The auxiliary 
condenser was constructed of two 9 . 5 mm (3/8 in) diameter, 
water-cooled copper lines helically coiled to a height of 
457 mm (18 in) . 

Cooling water for the test tube was provided, on a once- 
through basis, by two centrifugal pumps connected in series. 
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Figure 2.1. Schematic of Test Apparatus 
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Figure 2.2. Schematic Details of Test Section (Insert not Shown) 



Water to these pumps was provided through a large sump 

3 

(capacity about 0.4 m ), which was filled with a continuous 
supply of filtered tap water. The water could be throttled 
from zero to 0.55 1/s (8.8 gpm) , giving a maximum water 
velocity of 4.4 m/s (14.4 ft/sec) through the tube. 

A continuous supply of tap water was used for cooling 
the auxiliary condenser. Throttling the flow of tap water 
through the auxiliary condenser was the means used to vary 
the internal pressure of the test apparatus. 

A vacuum pump was used to create vacuum and for removal 
of noncondensing gases from the test section. The vacuum 
pump was in operation throughout the data collection time 
interval . 

B. SYSTEM INSTRUMENTATION 

The input power to the heaters was changed through a 
panel-mounted potentiometer. A converter with an input 
voltage of 440 VAC generated a signal, which was fed to the 
data acquisition system in order to calculate the power 
input to the heaters. A complete explanation of the con- 
verter function is given by Poole [Ref. 13]. 

The internal pressure of the system was measured manually 
by a U-tube, mercury-in-glass manometer graduated in milli- 
meters. Unavoidably, steam would condense in the manometer. 
Therefore, the varying height of the water column in the 
manometer needed to be accounted for when measuring the 
system pressure. 
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All temperatures were measured using 0.25 mm diameter, 
copper-constantan thermocouples. Two thermocouples were used 
for the steam, one for the condensate return and one for the 
ambient. The calibration procedure for these thermocouples 
is described in [Ref. 13]. The temperature rise of the cooling 
water through the test section was measured by two Hewlett- 
Packard (HP) 2804A quartz thermometers as well as by a 10- 
junction, series-connected, copper-constantan thermopile. 

The temperatures were fed directly into the data-acquisition 
system as described below. During this entire investigation, 
these two independent means of cooling-water-temperature rise 
measurements agreed to within ±0.02 K. Note that the coolant 
temperature rise was as low as 0.7 K for the smooth tube 
(under vacuum conditions), while it was in excess of 1.2 K 
for the finned tubes. Thus, the coolant temperature rise was 
measured with an uncertainty better than ±3 percent for finned 
tubes . 

C. DATA ACQUISITION SYSTEM 

An HP 3497A Data Acquisition System, as controlled by 
an HP 9826A computer, was used to monitor the system tempera- 
tures and the converter's signal representing the power input 
to the heaters . The data were immediately processed and also 
were stored on computer disks for later reduction and 
evaluation. 
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D. SYSTEM MODIFICATIONS 



For some of the early data taken near atmospheric pres- 
sure, the measured steam-side heat-transfer coefficient, 
when operating with a closed system, disagreed by as much as 
±30 percent between different data runs. This was later 
found to be due to the malfunctioning of the relief valve 
provided on the apparatus for safety reasons (see Figure 
2.1). This valve had a designed relief gauge pressure of 
52 mmHg. However, occasionally it was found to open with a 
system gauge pressure of only about 7 mmHg. Once opened, the 
valve remained open even if the system operating pressure 
was 10 mmHg below atmospheric, thus allowing air to enter 
the apparatus. To avoid this problem, subsequent runs were 
made without the relief valve and a new, more reliable, 
valve was ordered. In the meantime, extreme care had to be 
exercised to avoid system pressurization. 

Originally, the apparatus was provided with a 10.6 kW 
(3 ton) air-conditioning system to cool the cooling water 
sump. However, as also mentioned by Poole [Ref. 13], this 
system was inadequate to keep the sump temperature from 
rising. For this reason, this air-conditioning system was 
removed. By removing this unit, the cooling water friction 
losses were reduced and the maximum water velocity was in- 
creased from 4.2 to 4.4 m/s. In the previous stage of the 
investigation, fresh tap water was continuously fed to a sump 
while an equal amount of water was being drained, maintaining 
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a constant sump level. Since the sump was small, the hot 
water from the test section had to be returned to the sump. 

For this reason, there was a continuous variation of the 
sump water temperature. In order to solve this problem, the 
sump was replaced by an aluminum tank three times its size 
and the outlet hot water was drained directly after its 
exit from the test tube. 

The air ejector used to create vacuum in the test sec- 
tion could not provide a successful operation owing to the 
cycling of the feeding air pressure. Occasionally the pres- 
sure in the test section could become smaller than the pres- 
sure in the ejector's air line resulting in a back flow 
which could contaminate the apparatus. This problem was 
solved by the use of a vacuum pump. The suction line for 
this vacuum pump was a penetration in the auxiliary condenser. 
Unavoidably, there was a continuous removal of steam together 
with the non-condensing gases. To avoid steam from reaching 
the vacuum pump (vacuum pumps are designed to remove only 
air) , the mixture was passed through a heat exchanger where 
the steam was condensed. This heat exchanger was installed 
as shown in Figure 2.3. Note that no new supply of cooling 
water was used for this heat exchanger; instead, the tap water 
supply line into the water sump was used for this purpose. 

Careful and persistent investigations into the occurrence 
of partial dropwise condensation on the copper tube revealed 
that the presence of nylon components, such as tube bushings 
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Figure 2.3. Schematic Details of Purge System and Water Sump 



(see Figure 2.2) and numerous ferrules used in the apparatus, 
were a source of system contamination. This problem worsened 
when the system pressure was reduced for a vacuum run after 
the system had operated at near-atmospheric pressures. This 
was thought to be owing to the fact that outgassing increases 
with increasing temperatures and decreasing pressures. To 
solve these problems, new Teflon tube bushings were manufac- 
tured. Also, all nylon ferrules were replaced with Teflon 
ferrules. Note that Teflon has an outgassing rate about one- 
hundred times smaller than nylon [Ref. 14]. 

E. VACUUM INTEGRITY 

Vacuum tightness for any condensation experiment, 
especially under low pressure conditions, is of vital impor- 
tance. Any data taken while the apparatus is inadequately 
air-tight, will lead to false conclusions because of the 
deleterious effect of the non-condensing gases.' Many efforts 
were made in the past to achieve a leak-free system. During 
this thesis effort, there were leak, tests on a routine basis. 

A total of eight vacuum- tightness tests were carried out over 
a period of ten months, including the period devoted to taking 
data. The leak rate was always less than that represented 
by a pressure rise equal to 2 mmHg in a 24-hour period. 

These tests were carried out at a pressure of about 85 mmHg 
absolute . 
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F. TUBES TESTED 



An instrumented tube was fabricated from a thick-walled 
copper tube with an inner diameter of 12.7 mm (1/2 in) and 
an outer diameter of 18 mm (3/4 in) . The condensing length 
was 13.3 mm. This tube was provided with six wall thermo- 
couples, equally spaced around the tube periphery and 
located radially at the mid-point of the tube wall thickness. 
Details about the fabrication and instrumentation of this 
tube are given by Poole [Ref. 13]. Figure 2.4 shows a photo- 
graph of the instrumented tube. This tube was used to 
determine the leading constant of a Sieder-Tate type equation 
for the inside heat-transfer coefficient. 

A smooth tube of copper with dimensions similar to the 
instrumented tube was also fabricated. The purpose of the 
smooth tube was to determine the Sieder-Tate leading constant 
by using a modified Wilson plot, and to compare the data for 
the steam-side coefficient with the Nusselt prediction. 

Twenty one finned tubes were also fabricated. These 
tubes had the same inside diameter and length dimensions as 
the smooth tube. Further, they had a root diameter equal to 
the outside diameter of the smooth tube. The fins were 
rectangular in shape. The fin dimensions were systematically 
changed as shown in Table I. Figure 2.5 shows a photograph 
of the group of tubes with a fin thickness of 0.75 mm and a 
fin height of 1.0 mm. 
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Figure 2.4. Photograph of the Instrumented Tube 
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Photograph of a Set of Tested Tubes with 
Fins of t = 0.75 and e = 1.0 mm 





TABLE I 



Geometry of Finned Tubes 



Fin Spacing 
(mm) 




Fin 

Thickness 

(mm) 


Fin 

Height 

(mm) 


0.5, 1.0, 1.5, 


2.0, 


4.0 


1.0 


1.0 


0.5, 1.0, 1.5, 


2.0, 


4.0 


0.75 


1.0 


0.5, 1.0, 1.5, 


2.0, 


4.0 


0 . 5 


1.0 


1.0, 1.5, 


2.0, 


4.0 


1.0 


2.0 


1.0 






1.5 


1.0 


1.25 






1.25 


1.0 



G. INSERT 

Early investigations showed that with no internal en- 
hancement, the water-side thermal resistance could be as much 
as 90 percent of the overall resistance. Under such condi- 
tions, even a small discrepancy on the water-side resistance 
will translate into a substantially greater discrepancy on 
the condensing heat-transfer coefficient. Therefore, in 
order to improve the accuracy of the condensing heat- transfer 
coefficient, a decision was made to enhance the inside heat- 
transfer coefficient. 

For this purpose, a spiral insert was fabricated. The 
insert consisted of a 6.4 mm diameter stainless-steel rod 
together with a copper wire that was wrapped around the rod. 
The diameter and the pitch of the wrapped wire were 3.2 mm 
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and 20 mm respectively. After wrapping, the wire was properly 
soldered on the rod and then it was machined down to give a 
clearance of 0.5 mm between the outer wire diameter and the 
tube inside wall. This clearance was considered important 
to avoid heat conduction from the tube wall to the insert. 
Figure 2.6 is a photograph of the insert with one of the 
finned tubes. 
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Figure 2.6. Photograph of the Insert Together with 

a Finned Tube Tested 



Ill . 



SYSTEM OPERATION AND PROBLEMS 



Each tube to be tested was cleaned following the proce- 
dures listed in Appendix A. After installing a tube, the 
system was brought to operating pressure by following the 
procedures listed in Appendix B. The data were collected 
when steady-state conditions were achieved which were 
assumed to occur when the temperature of the steam in the 
test section and the inlet-outlet water temperatures were 
stabilized. Double data sets were taken for each flow 
setting starting with a cooling water flow rate of 80% 

(4.4 m/s). The next water flowmeter settings were 70, 60, 

45, 35, 26, and 20 percent followed by 55 and 80 percent. 
These settings were selected to provide nearly equal-spaced 
heat flux data points. 

Since the tubes were fabricated from copper, a highly 
conductive material, the time interval required to achieve 
steady-state conditions was not more than two minutes. A 
whole run (18 data sets) was completed in a time interval of 
less than one hour. After each change in the cooling water 
flow rate, the system pressure changed slightly. In order 
to keep this system pressure nearly constant, it was neces- 
sary to adjust the flow rate through the auxiliary condenser. 

During the early stages of this experiment, it was 
observed that the thermocouple readings were disturbed by an 
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external radio-frequency (RF) signal. This RF signal was 
generated by the special device used to control the power 
input to the boiler heaters. This device consisted of a 
silicon-controlled rectifier that chopped the input signal 
to obtain the necessary electrical power, while the chopped 
portion was emitted as an RF signal. This RF signal caused, 
for example, the steam temperature to fluctuate within ±4 K, 
while the thermopile readings fluctuated within ±0.5 K (note 
that the thermopile readings measuring temperature rise of 
cooling water can be as low as 0 . 5 K) . To alleviate these 
problems, all thermocouples were shielded against the RF 
signal using braided copper cables. This substantially 
reduced the temperature fluctuations. For example, after 
the shielding was installed, the thermopile readings fluctu- 
ated less than ±0.04 K. Further efforts revealed that the 
small fluctuations which persisted were due to the presence 
of a cable that carried an indication of the voltage to the 
boiler heaters. This voltage is of "DC" type and it was fed 
to the data acquisition system. Unfortunately, this voltage 
was not truly DC; instead it consisted of the RF disturbances 
described earlier. When this cable was disconnected, abso- 
lutely no fluctuations were observed. Thus, a switch was 
added on this cable to disconnect the source of the RF 
signal while thermocouple readings were being recorded. Soft- 
ware was modified to gather all thermocouple readings while 
this switch was turned off. When the computer was ready to 
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'read' the voltage to the boiler heaters, the program was 
coded to 'pause' and then remind the operator to turn on the 
switch. Immediately after taking this reading, the computer 
reminds the operator to turn off this switch again. 

The appearance of the condensate film was regularly 
observed through the viewport. In the event the condensate 
film appeared to be patchy, the data run was discontinued 
and the data already collected were discarded. During some 
runs, the film appeared to be totally filmwise, but the 
final data set indicated a condensing coefficient up to 10 
percent greater than the value for the initial data point. 
Recall that the initial and final data sets were for the 
same cooling water flow rate. It is worth noticing that it 
was easier to obtain filmwise condensation under vacuum 
conditions than at near-atmospheric pressures. This is 
probably due to lower vapor velocities associated with 
atmospheric runs (2 m/s for vacuum and 1 m/s for atmospheric 
conditions) . All the data presented in this thesis showed 
less than ±5 percent disagreement in condensing coefficient 
between initial and final data sets. 

A. THE DROPWISE CONDENSATION PROBLEM 

The main objective of this thesis effort was to take 
data under strict filmwise condensation conditions using 
copper tubes. Pure copper has low wettability with water in 
the presence of a minute contamination; thus, it has the 



47 



tendency to condense steam under the dropwise mode. Such a 
minute contamination is generally difficult to avoid. A 
slightly hydrophobic surface condition can lead to very 
false conclusions since the dropwise condensation mode of 
heat transfer is much more effective than the filmwise mode. 

In the first stage of this thesis effort, the tubes were 
prepared following the procedure recommended by Poole 
[Ref. 13]. This procedure included the cleaning of tubes 
with a mixture of sodium hydroxide and ethyl alcohol (equal 
amounts by mass) at a temperature of about 80 °C. This proce- 
dure was soon proven inadequate as the copper tube showed 
partial dropwise condensation within 20 minutes. This was 
verified during a four-hour endurance test. The water flow 
rate was held constant and every 30 minutes the temperature 
rise of the cooling water was recorded. This rise continu- 
ously increased and gradually the condensation changed to 
the dropwise mode. 

As a solution to this problem, it was decided to oxidize 
the tube surface more thoroughly using the above-mentioned 
mixture (50 percent) aqueous sodium hydroxide mixed with an 
equal part of ethyl alcohol) . The tube was painted with this 
mixture and was then heated with steam for about an hour 
with no water flowing inside the tube. This process resulted 
in the formation of a dark oxide layer on the copper surface, 
which had high wetting characteristics and an insignificant 
thermal resistance. After the oxidation, new endurance tests 
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were made. For a time period of three hours, the cooling 
water temperature rise was essentially unchanged, while the 
condensate film appearance was excellent. Steam-side coeffi- 
cients obtained with and without oxidation agreed to within 
±2 percent. This was probably due to the avoidance of even 
minute areas undergoing dropwise condensation. 

The most significant achievement of this oxidation was 
the reduction of data scattering (less than ±1 percent) , and 
the very good repeatability of similar runs (within ±3 percent) 

B. VAPOR VELOCITY LIMITATIONS 

The major assumption in the Nusselt theory is that the 
vapor velocity does not induce shear forces on the vapor- 
liquid interface. This assumption is not realistic in prac- 
tice. In order to regulate the pressure in the test section, 
it is necessary to have some steam condensed in the auxiliary 
condenser. This means that there must be steam flow into 
the auxiliary condenser. The steam flow rate (i.e., the 
vapor velocity) is a function of the power input to the 
heaters. In order to achieve low vapor velocities, the power 
input should be decreased. At low pressures, however, it 
was evident that the fluctuation of system pressure increased 
with a decrease in boiler power. This fluctuation was due 
to the intermittent breakup of vapor bubbles on the water- 
vapor interface in the boiler. Surprisingly, this small 
fluctuation in system pressure, and hence vapor velocity. 
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reflected in the cooling water outlet temperature. To 
minimize this fluctuation, a vapor velocity of approximately 
2 m/s past the tube was selected for vacuum runs. For the 
atmospheric runs (maximum power input) , the cooling water 
outlet temperature fluctuation was less than ±0.005 K, while 
for the vacuum runs (half of maximum power input) , the fluc- 
tuation was less than ±0.05 K. 

A low vapor velocity is also undesirable from the film- 
wise point of view. When the vapor flow is low, then the 
condensation rate decreases and the film can very easily 
break, resulting in a change of the condensation mode. Once 
the film has broken, there is no possibility for the tube to 
resume the filmwise mode. Also, when the vapor velocity is 
high, there is a very strong wash-off effect on the tube, 
resulting in the reduction in the tendency of the tube to 
change the condensation mode. 

It was decided to use 2 m/s and 1 m/s vapor velocities 
for the vacuum and atmospheric runs respectively. These 
velocities are considered much smaller than the actual ones 
in condensers. However, due to the limited power input, the 
1.0 m/s used for the atmospheric runs was indeed the maximum 
possible. For the vacuum runs, the 2.0 m/s vapor velocity 
was achieved by using half of the available power. Some 
attempts were made to increase the vapor velocity to 4.0 
m/s, but it was considered unsafe due to the violent and 
rapid explosions of the vapor bubbles in the boiler. 
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IV. DATA REDUCTION 



During this thesis effort, all the programs, including 
property functions and the calibration curves for the 
cooling water flowmeter and for all thermocouples, were 
carefully checked. Only a very few minor corrections had to 
be made. The programs have the following noteworthy features 

1. In order to obtain a meaningful mixing-cup temperature 
at the coolant outlet, a mixing chamber was provided at the 
outlet as shown in Figure 2.2. As also discussed by Poole 
[Ref. 13], the considerable pressure drop created by this 
mixing chamber introduced a temperature rise in the cooling 
water. With no steam condensation taking place, the coolant 
temperature rise was measured for various water velocities. 
The results were plotted as shown in Figure 4.1. When taking 
condensation data, the coolant temperature rise was corrected 
by subtracting the temperature rise due to this frictional 
heating at the mixing chamber. For this purpose, the least- 
squares line shown in Figure 4.1 was used. 

2. The smooth end sections of the tube, which were neces- 
sary to hold the tube in place (see Figure 2.2), were acting 
as axial fins. Since copper has a very large thermal 
conductivity, axial heat conduction occurred. This heat 
conduction was incorporated into the data reduction process 
by computing the fin efficiency <bf the smooth end sections 
of the tube. 
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Figure 4.1. Mixing Chamber Calibration Curve 



3. The tap water temperature changes considerably depend- 
on the season: it is about 15 °C in winter and about 

22 °C in summer. Because of these changes, it was necessary 
to incorporate water property variations in calibrating the 
flowmeter. For this purpose, the following function was 
used, as done so in [Ref. 15] , to calculate a correction 
factor for viscosity changes. 



C f = 1.0365-1.9644(10 3 ) Tin + 5 . 25 ( 10 6 )Tin 2 (4.1) 



The value of the above function at the time of flowmeter 
calibration (i.e., at Tin = 22.0 °C) was 0.995. Thus, in 
order to obtain the actual mass flow rate for an inlet 
temperature of Tin, the following equation was used: 



m 



act 



= m 



calc 0.995 



(4.2) 



The major steps for the data-reduction programs are 
listed below. 

A. PROGRAM SIEDER 

This program calculates the constants of the Sieder-Tate 
type equation from the instrumented-tube data. 

1. Read the cooling water inlet and outlet temperatures 
and compute the average temperature. 

2. Read the six wall thermocouples and compute the 
average wall temperature. 
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3. Convert the water volume flowmeter setting to mass 
flow rate by dividing the volumetric water flow rate with 
the water velocity. 

4. From continuity, compute the cooling water velocity. 

5. Correct the outlet water temperature for the mixing 
chamber effect and then re-compute the average temperature. 

6. Compute the amount of heat transferred to the cooling 
water . 

7. Compute the wall thermal resistance and then find the 
inner wall temperature (from step 2, the average temperature 
at the radial location of the thermocouples, is known) . 

8. Compute the log-mean- temperature difference. 

9. Compute the efficiency of the left and right tube 
extensions . 

10. Compute the inside heat- transfer coefficient. 

11. Compute the Nusselt number. 

12. Compute the leading constant of the Sieder-Tate 
coefficient for this data set. 

13. Repeat steps 1 to 12 for all the data sets of the run. 

14. Make a least-squares straight line fitting of the 
form: 



ao i/-> P 0.14 , . 
r-r _ 0.8 „ 1/3 , c, + B 

Nu = C Re Pr ( — ) 



(4.3) 



w 



The C is the unknown leading constant, while B is an addi- 
tional constant found to give an improved fitting. 
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B. MODIFIED WILSON PLOT PROGRAM 



This program calculates the leading constant of the 
Sieder-Tate equation from the smooth-tube data. 

1. Assume a value for the leading constant of the Sieder- 
Tate equation. 

2. Make the steps 1, 3, 4 and 5 of previous program 
"SIEDER. " 

3. Compute the heat flux from the specific heat of water, 
the mass flow rate, the cooling water temperature difference 
and the outside tube surface. 

4. Compute the log-mean- temperature difference between 
the steam temperature and the inlet and outlet cooling water 
temperatures . 

5. Compute the overall heat-transfer coefficient. 

6. Assume an outer tube surface temperature in order to 
determine the film temperature for the evaluation of conden- 
sate properties. 

7. Calculate the outside heat- transfer coefficient by 
using the Nusselt equation. 

8. Compute the outer tube surface temperature. 

9. Iterate between steps 7 and 8 until two successive 
outer surface wall temperatures are within ±1 percent. 

10. Assume a viscosity correction factor for the 
Sieder-Tate equation. 

11. Compute the fin efficiency of the extending parts of 
the tube. 
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12. Compute the inside heat- transfer coefficient. 

13. Compute the inner surface temperature in order to find 
the film temperature on the coolant side. 

14. Compute the viscosity correction factor and iterate 
between steps 11 through 14 until two successive viscosity 
correction factors are within ±1 percent. 

15. Compute the new Sieder-Tate coefficient and iterate 
between steps 2 through 15 until two successive Sieder-Tate 
coefficients are within ±0.5 percent. 

16. Compute the Sieder-Tate coefficient for all the data 
sets of the run. 

17. Make a least square straight line fitting to establish 
the final leading constant that satisfies all the data sets. 

C. MAIN DATA REDUCTION PROGRAM ( DRP4 ) 

This program was used for computing the condensing heat- 
transfer coefficient for both finned and smooth tubes. 

1. Make steps 1 through 5 of program "SIEDER. " 

2. Compute the log-mean temperature difference between 
the steam temperature, and the inlet and outlet cooling 
water temperatures . 

3. Compute the heat flux. 

4. Compute the overall heat-transfer coefficient. 

5. Compute the wall conduction resistance. 

6. Compute the inside heat- transfer coefficient. 

7. Compute the condensing heat-transf er coefficient by 
subtracting the wall and the inside resistances from the 
overall resistance. 
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V. RESULTS AND DISCUSSION 



Numerous data runs were made using the procedures 
described in Chapter III. Each tube was tested at least 
three times, on different days, in order to obtain repeata- 
ble results. Data were taken for all the tubes at both 
near-atmospheric and vacuum (85 mmHg , 1.7 psia) conditions. 
Complete filmwise condensation was maintained for the data 
runs and the mass concentration of non-condensing gases was 
held between 0.0 and -1.5 percent during all the testings. 
The negative value was indicative of either a slight super- 
heat or an inaccurate manometer reading. It was estimated 
that an error of ±0.2 K in steam temperature or an error of 
1.0 mmHg in system pressure translates into approximately 
±1.2 percent error in the non-condensing gas concentration. 
Thus, the non-condensing gas concentration was essentially 
zero to within the accuracy of the measurements. 

A. INSIDE HEAT-TRANSFER COEFFICIENT 
1 . Instrumented Tube Results 

The instrumented tube was used to determine the 
inside heat-transf er coefficient mainly owning to the use of 
the insert to enhance the heat- transfer performance. By 
the use of the instrumented tube, the wall temperature was 
determined after averaging the six thermocouple readings. 
Therefore, the main assumption in the calculation of the 
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inside heat-transf er coefficient is that the wall tempera- 
ture profile is uniform. However, Figures 5.1 and 5.2 show 
the experimentally-measured temperature profiles around the 
tube, for the maximum and minimum cooling water flow rate 
settings, for the vacuum and near-atmospheric runs, respec- 
tively. It is obvious that the wall temperature profile is 
not uniform; therefore, this direct averaging procedure may 
introduce an uncertainty the magnitude of which is unknown. 
However, owing to the time constraint, further evaluation 
of this uncertainty was considered to be beyond the scope 
of this investigation. These figures also show least-squares- 
fit cosine curves of the form; 

~ = 1 - a cos - (5.1) 

AT 71 

It is worth noticing that the temperature profile corres- 
ponding to the high heat flux is flatter than the profile 
corresponding to the low heat flux. Also, the temperature 
profile of the near-atmospheric run is flatter than the 
temperature profile of the vacuum run. It is concluded that 
as the heat flux increases, hi increases, hence ho/hi 
decreases and the temperature profile becomes more uniform; 
therefore, the assumption of an isothermal wall becomes more 
reasonable. This is one more reason why the internal 
enhancement (use of insert) improves the accuracy of the 
data . 
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Figure 5.1. Temperature Distribution for Instrumented 
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Figure 5 . 2 . Temperature Distribution for Instrumented 



Tables II and III give the values of the amplitude of the 
cosine fit for all data points taken, under both pressure 
conditions . 



TABLE II 

Instrumented Tube Results for the Vaccum Run 



Heat Flux 
(MW/m 2 ) 


ho 

(kW/m 2 .K) 


ho/hi 


a 


0 . 209 


10.89 


0 . 301 


0.135 


0 . 207 


11.28 


0.348 


0.142 


0.199 


11.45 


0.430 


0.158 


0.191 


11.61 


0.507 


0.168 


0 .184 


11.76 


0.567 


0.178 


0.173 


12.05 


0.692 


0.185 


0.160 


12.09 


0 . 832 


0.202 





TABLE 


III 






Instrumented 


Tube Results 


for 


the Atmospheric Run 


Heat Flux 
( MW/m 2 ) 


ho 

(kW/m 2 .K) 




ho/hi 


a 


0.488 


9.18 




0 . 201 


0.115 


0.486 


9.19 




0 . 221 


0.119 


0 . 479 


9.25 




0 . 249 


0.133 


0.462 


9.38 




0 .289 


0.137 


0.447 


9.49 




0.346 


0.149 


0.417 


9.75 




0.437 


0.164 


0.384 


10.05 




0.585 


0 . 179 
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Also the same tables contain additional information 



about the heat flux, the outside heat- trans fer coefficient 
and the ratio of outside to inside heat-transfer coefficients 
for all the cooling water flow rate settings. 

By using the average reading of all thermocouples 
for the wall temperature, the complete form of the Sieder- 
Tate-type equation was determined to be : 

Nu = 0.0635 Re 0 * 8 Pr 1/3 (-) 0 ' 14 + 26.4 (5.2) 

^w 

The constant B = 26.4 is an additional constant found to 
improve the fitting. Note that the leading coefficient is 
about 2.5 times greater than the well-known Sieder-Tate 
coefficient of 0.027 [Ref. 16] for plain tubes. This signi- 
ficant difference is mainly due to the coiled insert used 
to boost the inside heat- transfer coefficient. A discussion 
of the advantages and disadvantages of the coiled insert was 
given in Chapter II, Section G. Figure 5.3 shows the data 
points for both pressure conditions together with a fit of 
equation (5.2) . 

It is worth mentioning that the use of thermocouples 
to determine the wall temperature may introduce yet another 
uncertainty (in addition to the incorrect averaging mentioned 
earlier) in the measurement of each local wall temperature. 
This is due to the completely different thermal conductivity 
of the wall area in the immediate vicinity of the thermocouple 
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Figure 5.3. Correlation of the Water-Side Coefficient 

Obtained with Instrumented Smooth Tube 



hot junction. For example, the thermocouple wires are 
covered with plastic; therefore, they present a low thermal 
conductivity. Also, in order to embed the thermocouple, 
there must be sufficient clearance, which certainly presents 
an additional thermal resistance to the heat flow. 

The advantages of the use of the instrumented tube 

are : 

1 . It does not depend on other theoretical formulations 
that may have unavoidable errors owing to the assumptions 
made to formulate the problem (for example, the Modified 
Wilson Plot technique uses the Nusselt-type equation for the 
outside heat-transfer coefficient) . 

2. It is not too sensitive to vapor shear effects since 
these effects are lumped into the thermocouple readings 
(i.e., the vapor shear reduces the thickness of the conden- 
sate layer; therefore, the thermocouples would give higher 
temperatures) . 

A computer program, named "SIEDER, " was written to 
reduce the data taken on the instrumented tube. For com- 
pleteness, the raw data taken on the instrumented tube for 
both pressure conditions are given in Appendix C, Table VI. 
Hereinafter, the data reduction method that uses the instru- 
mented tube water-side heat-transfer coefficient will be 
referred to as the "Direct" method. 

2 . Uninstrumented Tube Results 

The uninstrumented smooth-tube data were used to 
determine the inside heat- transfer coefficient by the use of 
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the modified Wilson plot technique. Hereinafter, this tech- 
nique will be referred to as the "Wilson" method. Also the 
smooth tube results concerning the outside heat-transfer 
coefficient were used as a basis for comparison with the 
finned- tube results. 

The analysis procedures used for the Wilson method 
is the same as that given in [Ref. 17] . While a Sieder- 
Tate-type equation was used to represent the inside heat- 
transfer coefficient, a Nusselt type expression was used for 
the outside coefficient as given by: 



h 



o 



a[ 



k fpf(p f p v )h fg . 1/3 



y f D o q 



■3 



(5.3) 



where a is a disposable constant that must be found itera- 
tively. Note that with zero vapor velocity, the Nusselt 
prediction results in a value of 0.655 for a. However, 
since it is impossible to obtain zero vapor velocity, a 
higher value is generally obtained. During this work, an 
approximate value of 0.72 was obtained for a. 

A computer program, named "WILSON," was written to 
run the Wilson method on the smooth-tube data. A total of 
four runs (two each under vacuum and at near-atmospheric 
pressures) were made. These runs resulted in values of the 
Sieder-Tate coefficient of 0.0716, 0.0714, 0.0700 and 
0.0708, with an average value of 0.071. This value is about 
10 percent higher than the value found with the instrumented 
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tube; i.e., the Wilson method results in a 10 percent 
smaller water-side resistance (less conservative) than the 
resistance predicted by the use of the instrumented tube. 
Again for completeness, the raw data taken on the uninstru- 
mented tube under both pressure conditions are given in 
Appendix C, Table V. 

3 . Summary of the Inside Heat-Transfer Coefficient 

In order to compare the two above-mentioned methods 
(Direct and Wilson) , the data taken on a smooth tube under 
vacuum conditions were reduced to yield the outside heat- 
transfer coefficients using the water-side coefficients 
derived by each method. For this purpose, these data were 
plotted in Figure 5.4 together with the original Nusselt 
equation for a single horizontal tube, re-written in terms 
of the indicated non-dimensionalized terms: 

Nu l /q 

— = 0.725 F /r> (5.4) 

Re 

Also shown in this figure is the correlation of Fujii, et 
al [Ref. 18] : 



, n % = 0.96 F 1/5 (5.5) 

Re 0 ' 5 

The Fujii et al . correlation is based on numerous data 
taken on tube bundles in both in-line and staggered arrange- 
ments . It can be seen that the coefficient on the right-hand 
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Figure 5.4. Comparison of Smooth-Tube Data with Fujii Correlation 



side of equation (5.5) is larger than the value in equation 
(5.4) owing to the presence of vapor shear (recall that the 
Nusselt equation is derived for zero-vapor-shear conditions) . 
For completeness, the error bands of the data used for this 
correlation are also shown. 

Note in Figure 5.4 that the data of this investiga- 
tion on the smooth tube are slightly higher than those 
predicted by equation (5.5) . Also, it is worth noting that 
the use of the Wilson method results in data that are in 
better agreement with the Fujii correlation than the data 
given by the use of the Direct method. Thus, with an appro- 
priate function for the outside heat-transfer coefficient 
and accurate smooth- tube data, the modified Wilson plot 
technique appears to offer a satisfactory and simple method 
for obtaining the inside heat- transfer correlation. However, 
since the error band of the Fujii correlation is much wider 
than the disagreement between this correlation and the 
present data, it may be unwise to draw any conclusions as 
to the accuracy of these two methods (Direct or Wilson) . 

As a result, it is not possible to say which method of measur- 
ing the inside heat-transfer coefficient is more accurate. 
However, except where noted below, all the outside heat- 
transfer coefficients obtained in this thesis were calcu- 
lated using the Direct method. 
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B. PERFORMANCE OF FINNED TUBES 



Data were taken for the finned tubes given by Table I, 
under vacuum (85 mmHg) and at near-atmospheric pressures 
with vapor velocities of 2 m/s and 1 m/s respectively. The 
computed steam-side coefficient was based on the smooth tube 
outside area (i.e., the area if the fins were machined off) . 

The data reductions were performed using the program 
named "DRP4." This is a revised version of the program 
(DRP2) used by Poole [Ref. 13]. The raw data taken on finned 
tubes are presented in Appendix D. 

1 . Condensation on Fins in General 

Figure 5.5 shows the variation of heat flux with the 
steam-side temperature drop for a set of five finned tubes 
under vacuum conditions. Both the thickness and height of 
the fins were 1.0 mm, while the spacings were 0.5, 1.0, 1.5, 
2.0 and 4.0 mm. Notice that the slope of the curves in- 
creases with fin spacing up to the point where the optimum 
spacing (1.5 mm) is reached and then the slope decreases 
with further increase in fin spacing. This means that the 
enhancement of the condensation process is highly dependent 
on the fin spacing. The data were fitted with least-squares 
curves using the following relation: 

q = a AT b , (5.6) 

and AT was computed by using the following relation: 
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Figure 5.5. Variation of Heat Flux with Steam-Side Temperature Drop 



q 



h AT 
o 



(5.7) 



where q was measured experimentally, while h Q was obtained 
by subtracting inside and wall resistances from the measured 
overall resistance. From equation (5.6), it is possible to 
calculate the heat flux for any temperature difference be- 
tween the steam and the tube wall when the constants a and 
b are known . 

Table IV gives the values of a and b for all the 
tubes tested for both pressure conditions. Notice that the 
value of the constant a has almost the same range for all 
the finned tested tubes for both pressure conditions. The 
values of both constants a and b increase as the fin spacing 
increases up to the spacing of 1 . 5 mm (it will be shown 
later that this is the optimum spacing) and then they 
decrease. The range of the exponent b is different between 
the vacuum and near-atmospheric runs; i.e., the range for 
the vacuum runs is 0.76 to 0.89, while for the near-atmospheric 
runs it is 0.87 to 0.94. 

Figure 5.6 shows a relative plot of the outside 
condensing coefficient for the same set of tubes just 
described. For comparison purposes, the smooth- tube data 
and the Nusselt prediction are plotted on the same figure. 

All tubes tested show good enhancement over the smooth tube. 

It is worth noticing that even the smooth tube results are 
higher than the Nusselt prediction; this is mainly 
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with Heat Flux 



attributable to Nusselt's conservative assumption which 
neglects vapor shear. 

Combining equations (5.6) and (5.7) , h Q can be 
written as: 



1/b (b-l)/b 

o - a q 



(5.8) 



It appears that equation (5.8) results in curves almost 
parallel to the Nusselt prediction, indicating that the 
effect of heat flux is the same on the finned surface as on 
the smooth tube . 

On the above figure, the dotted lines represent the 
fittings of the data when the data were processed by the 
Wilson method. Notice that the condensing coefficients 
given by the Wilson method are lower by about 10 percent 
than those given by the Direct method (recall that the 
water-side resistance given by the Direct method is higher 
than that given by the Wilson method) . 

It is also shown in the above figure that as the 
heat flux decreases, the vapor-side heat-transfer coeffi- 
cient increases. Note that there is about a 10 percent 
difference between the high and the low heat flux results. 

The reason for this increase is that, as the flux decreases, 
the condensate thickness decreases, resulting in less 
conduction resistance and improved heat-transfer coefficient. 
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2 . Significance of Fin Spacing 

This section presents results showing the variation 
of the condensing coefficient with heat flux, as well as the 
variation of heat flux with steam-side temperature drop for 
both pressure conditions, for the four sets of finned tubes. 

Figures 5.10 to 5.25 show that for all the fin thicknesses 

% 

(i.e., 1.0, 0.75 and 0.5 mm) and for each of the fin heights 
(i.e., 1.0 and 2.0 mm), the best performance was obtained 
with the tubes of 1.5 mm fin spacing. This was true for 
both pressure conditions. The second best performance was 
given by the tubes with spacings of either 1.0 mm or 2.0 mm. 
The third best performance was given by the tubes with fin 
spacing of 4.0 mm and the tubes with a fin spacing of 0.5 mm 
showed the poorest performance. The performances of the 
last two tubes however were interchanged as the pressure 
conditions were changed from vacuum to near-atmosph eric . An 
effort to explain this behavior will be made in the next 
sections, where the influence of the fin thickness and 
height on the heat-transfer performance is discussed. 

As shown in Appendix E (Error Analysis) , the error 
in the calculation of the condensing coefficient increases 
as the heat flux decreases. For this reason, the comparison 
of finned tubes should be done in the regime where the error 
is the minimum possible. This means that the data points 
corresponding to the high-heat-flux settings should be 
chosen. The comparison of the finned tubes is made through 
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the enhancement ratio, which is defined to be the ratio of 
the outside heat-trans fer coefficient of a finned tube to 
the smooth-tube value (both evaluated on the same heat 
flux) . However, the capability of the smooth tube to remove 
heat is limited (see Figure 5.6 showing that the maximum 
heat flux on the smooth tube is only about 30 percent of the 
maximum heat flux on the best finned tube) . Thus, it was 
necessary to extrapolate the smooth-tube data to higher heat 
fluxes in order to make the comparison of the finned tubes 
on an equal heat flux basis. 

2 

Choosing 0.35 and 1.0 MW/m as reference heat fluxes 
for the vacuum and near-atmospheric runs, respectively, 
cross plots of the enhancement ratio versus the fin spacing 
are given in Figure 5.7 for the set of tubes with 1.0 mm fin 
thickness and height. It is shown that the optimum fin 
spacing found is 1.5 mm. The 0 . 5-mm fin spacing which pro- 
vides the most area increment has one of the worst performances. 
This poor performance can be easily explained by the fact 
that the entire interfin area was covered by the condensate. 
Visual observations showed that the condensate retention 
angle was 180 degrees, which is also very well predicted by 
equation (1.12) . The thermal resistance induced by this 
thick condensate layer deteriorates any benefit from in- 
creased area; therefore, the heat-transfer performance is 
highly reduced. 
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Figure 5.7. Cross Plot Showing Best Fin Spacing Among 

Tubes with t - 1 mm and e = 1 mm 



As the fin spacing increased from 0.5 to 1.0 mm, 
the retention angle decreased and a miniscus in the interfin 
area at the top of the tube appeared (see Figure 5.8(a)) . 

For the 1.0 mm fin spacing, the retention angle was about 
110 degrees. The change of the flooding characteristics is 
the reason for the sudden increase of the enhancement ratio 
as the fin spacing changed from 0.5 to 1 mm. Visual obser- 
vations showed that the optimum fin spacing (1.5 mm) resulted 
in a condensate film profile (in the interfin area at the 
top part of the tube) similar to that shown in Figure 5.8(b) . 
Notice that the bottom of the miniscus is just touching the 
tube surface. When the fin spacing is further increased, 
the miniscus radius becomes larger and the retention angle 
reduces even more. At the 2.0 mm fin spacing, the miniscus 
is no longer an arc of a circle (see Figure 5.8(c)). 

As mentioned in the beginning of this chapter, the 
condensing-coef f icient calculation is based on the smooth 
tube external diameter; therefore, the effect of area addi- 
tion is lumped together with all the other factors that may 
influence the heat- transfer process. In order to give some 
insight to the effect of all the other factors, the enhance- 
ment ratio was divided by the area gain (i.e., the ratio of 
the condensing area of finned tube to the area of smooth 
tube) of each finned tube in order to obtain the normalized 
enhancement ratio. Figure 5.9 shows this normalized enhance- 
ment ratio. It is clear that the enhancement is due to 
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(a) 




(b) 




(c) 



Figure 5.8. 



Miniscus Profile for Various Fin Spacings 
(a) s = 1.0 mm, (b) s = 1.5 mm, (c) s = 2.0 mm 
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more than just an area gain. This verifies the measure- 
ments of Yau et al. [Ref. 3], and an explanation for this 
behavior was presented in Chapter I. Also, Figure 5.9 shows 
that the normalized enhancement ratio increases as the 
spacing increases and that the rate of this effect gets 
smaller for fin spacings larger than 1.5 mm. This reduction 
is attributable to the fact that the pulling action (thinning 
effect) diminishes as the fins get further apart. Equation 
(5.9) gives that the pressure gradient between points C and 
B (see Figure 1.1) is inversely proportional to the increase 
of fin spacing (s). 



AP 

dP \ ^ CB 

dx I s/2 

' CB 



(5.9) 



Therefore, the thinning effect decreases as the fin spacing 
increases. It. can be concluded that any assumption that 
neglects the surf ace- tens ion effect may not be realistic 
when modelling the steam condensation process on fins . 

Careful observation of Figure 5.7 shows that the 
performance of the tubes of 0 . 5 and 4 . 0-mm fin spacing 
interchanges between the vacuum and atmospheric runs, 
respectively. The heat that was removed during the vacuum 
runs was about three times smaller than the heat removed 
during the near-atmospheric runs; thus, the rate of conden- 
sate generation for the vacuum runs was three times smaller 
than for atmospheric runs. For all the vacuum runs, no 
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matter which tube was used, the condensate flow pattern was 
invisible. During the near-atmospheric runs, the flow 
pattern was also invisible except for the tubes of 0.5 mm 
fin spacing. In these cases (0.5 mm fin spacing), the 
draining rate was very fast and ripples in the interfin 
space were clearly visible. If this is the case, then the 
convective action of the condensate layer should enhance 
the condensation process. After this observation, it became 
clear why the performance of the tubes with 0.5 mm spacing 
became better than the performance of the tubes with the 
4.0 mm spacing, for the near-atmospheric runs. 

As mentioned earlier, Honda and Nozu [Ref. 9] 
present a detailed comparison of various theoretical models 
to predict the condensing coefficient on finned tubes. All 
the models have the assumption that the inertia term in the 
momentum equation and the convective term in the energy 
equation can be neglected in comparison with the other 
terms. The author feels that the above assumption is 
realistic for all the fin spacings except those that corres- 
pond to completely flooded tubes. To further support this 
argument, the reader is referred to [Ref. 9] where comparison 
of the four models is presented with various experimental 
data, among which are the steam data given by Rose [Ref. 3]. 
Honda and Nozu show that the maximum disagreement between 
the data of Rose and predictions (all four models) is for 
the fin spacing of 0.5 mm (completely flooded tube) . 
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Figure 5.10. Variation of Condensing Coefficient with Heat 

Flux for the Set of Tubes with t = 1.0 and 
e = 1.0 mm (Vac. Runs) 
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Figure 5.12. Variation of Condensing Coefficient with Heat 

Flux for the Set of Tubes with t = 1.0 
and e = 1.0 mm (Atm. runs) 
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Figure 5.13. Variation of Heat Flux with Steam-Side Temperature 

Drop for the Set of Tubes with t = 1 and e = 1 mm 
(Atm. Runs) 



2.0 m/s 
85 mmHg 




O Q 

(0 ID 



□ □ Q 

^ CO C\l 

OH / °M 



□ 



O 



□ 



87 



Figure 5.14. Variation of Condensing Coefficient with Heat Flux 

for the Set of Tubes with t = 0.75, e = 1.0 mm 
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Figure 5.15. Variation of Heat Flux with Steam-Side Temperature 

Drop for the Set of Tubes with t = 0.75 and e = 1.00 
mm (Vac. Runs) 



09 







ID ^ CO CM 

OH / 0 H 



r\ 

CM 

E 

\ 



\ 

0" 



<D 

5 



89 



Figure 5.16. Variation of Condensing Coefficient with Heat Flux for 
Set of Tubes with t = 0.75 and e = 1.0 mm (Atm. Pirns) 
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Figure 5.17. Variation of Heat Flux with Steam-Side Temperature 

Drop for the Set of Tubes with t = 0.75 and 
e = 1.0 mm (Atm. Runs) 
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Figure 5.18. Variation of Condensing Coefficient with Heat Flux 

for the Set of Tubes with t = 0.5 and e = 1.0 mm 
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Figure 5.19. Variation of Heat Flux with Steam-side Temperature 

Drop for the Set of Tubes with t = 0.5 and e = 1.0 
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Figure 5.20. Variation of Condensing Coefficient with Heat 

Flux for the Set of Tubes with t = 0.5 and 
e = 1 . 0 mm (Atm. Runs) 
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Figure 5.21. Variation of Heat Flux with Steam-Side 

Drop for the Set of Tubes with t = 0.25 
e = 1.0 mm (Atm. Runs) 
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Figure 5.22. Variation of Condensing Coefficient with Heat Flux 

for the Set of Tubes with t = 1.0 and e = 2.0 mm 
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Figure 5.23. Variation of Heat Flux with Temperature Drop 

for the Set of Tubes with t = 1.0 and e = 2.0* mm 
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Figure 5.24. Variation of Condensing Coefficient with Heat Flux 

for the Set of Tubes with t = 1.0 and e = 2 . 0 mm 
(Atm. Runs) 
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Figure 5.25. Variation of Heat Flux with Steam-Side Temperature 

Drop for the Set of Tubes with t = 1.0 mm and 
e = 2 . 0 mm (Atm. Runs) 



3 . Significance of Fin Thickness 



Three sets of finned tubes were manufactured with 
fin thicknesses of 0.5, 0.75 and 1.0 mm, while the fin 
height was held constant at 1.0 mm. Each set consisted of 
five tubes with fin spacings of 0.5, 1.0, 1.5, 2.0 and 4.0 

mm. The relative plots of the outside heat-transfer coeffi- 

* 

cient for each group of tubes, tested under vacuum and 
atmospheric conditions, were already shown in Figures 5.10 
to 5.21. Notice that the shape of the curves to fit the 
data points is very nearly independent of the fin thickness. 
Figure 5.26 is a cross plot showing the enhancement ratio as 
a function of the fin spacing for both pressure conditions, 
with the fin thickness as a parameter. Notice that for 
almost all the fin spacings, the 0 . 5-mm fin thickness gave 
the smallest gain even though it provides the maximum 
condensing area. An effort to explain this phenomenon will 
be made at the end of this section. Between the 0.75 and 
1.0-mm fin thicknesses, it is not very clear which of the 
two is superior. For the vacuum runs with fin spacing 
around the optimum value (i.e., 0.5 mm < s < 1.5 mm) it 
appears that the 1.0-mm thickness is a stronger candidate 
for the optimum fin thickness. On the other hand, for 
atmospheric runs with most fin spacings, the 0.75-mm fin 
thickness appears to be the best candidate. 

In order to further study the dependence of fin 
thickness on the enhancement ratio, data shown in Figure 
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Figure 5.26. Cross Plots Showing Best Fin Spacing Among Sets 

of Tubes with t = 0.5, 0.75 and 1 mm 



5.26 were cross plotted as shown in Figure 5.27. As mentioned 
earlier, initially there were only three tubes with fin 
thickness of 0.5, 0.75 and 1.0 mm. For vacuum conditions, 
the inadequacy of this fin thickness range can be seen 
since the optimum thickness has not been bracketed. There- 
fore, an additional tube with a fin thickness of 1.5 mm, 
a fin height of 1.0 mm and a fin spacing of 1.0 mm, was 
manufactured. The enhancement ratio found on this tube is 
also plotted in Figure 5.27. It can be seen that the tube 
with a 1.0-mm fin thickness is still the best candidate for 
the vacuum runs. This figure also shows possible error 
bands for the computed enhancement ratios. Since the 
possible error bands are in excess of the enhancement ratios 
at these two fin thicknesses, no attempts were made to 
explain why the optimum fin thickness is affected by the 
pressure condition. Therefore, it can be concluded that the 
optimum fin thickness is around 0.75 to 1.0 mm for both 
pressure conditions. 

As the fin thickness increases with fin spacing held 

constant, the fin density (fins/mm) decreases. In the 

limit, the tube will become a smooth tube with an outer 

diameter qual to the fin diameter. Recall that in the 

Nusselt equation (5.3), the tube diameter is in the denominator; 

therefore, the condensing coefficient is inversely propor- 

1/4 

tional to the D q for constant wall temperature, or, at 

a constant wall heat flux, the condensing coefficient is 

1/3 

inversely proportional to D ' . Thus, 
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Figure 5.27. Influence of Fin Thickness on Enhancement Ratio 

for Tubes with p = 1.0 mm, for both Pressure Conditions 



(5.10) 
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1/3 



where h' is based on the outer diameter of the tube, 
t-*» 

including the fin height. Now, if this condensing coeffi- 
cient is based on root diameter (as it is done throughout 
this thesis) , then 
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(5.11) 



The combination of equations (5.10) and (5.11) results in: 



D 2/3 

h , = h (~) 

t-*» o D 

o 



(5.12) 



This means that there will be very little enhancement by 

increasing the fin thickness to the limit. For example, 

with D- = 21.0 mm and D =19.0 mm, the ratio bn /h 
f o t+°° o 

would equal 1.07. 

When the fin thickness is reduced, there is a finite 
number of fins that can be added in a given tube length 
(i.e., it is not physically possible to take a smooth tube 
when fins of infinitesimal thickness are present) . Since 
there is a significant increase in heat- transfer area as the 
fins become thinner, an enhancement for the condensing 
process may be expected on the basis of area increment. 
Unfortunately, this was not the case as discussed below. 
Equation (1.12) predicts the retention angle and it is shown 
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that the condensate retention is independent of fin thick- 
ness. This was visually verified during this investigation. 
Comparing two tubes which have the same length and fin 
spacing, the one that has thinner fins would hold more 
condensate when the retention angle is the same. This means 
that part of the metal of a thick fin is replaced by conden- 
sate, or in other words, the thin-finned tube presents a 
smaller effective thermal conductivity than the thick-finned 
tube resulting in lower performance. Also, as the fin 
thickness decreases, the fin efficiency may play a role for 
tubes made of low- thermal-conductivity metal. However, for 
copper tubes with fins of 0.5 mm thickness and 1.0 mm 

height, the fin efficiency was estimated as 99 percent, 

2 

when the outside heat-transfer coefficient was 40 kW/m • K. 

Since both fin thickness extremes show that the 
heat-transfer process is not enhanced, it appears that the 
optimum thickness found (around 0.75 to 1.0 mm) is reasonable. 

Careful observation of Figure 5.26 reveals that the 
tubes with 0 . 5-mm fin spacing do not follow the trend of the 
tubes with larger fin spacings. In this case, the maximum 
performance was given by the tube with the 0.75-mm thick- 
ness, then by the tube with 0 . 5-mm thickness and finally by 
the tube with 1.0-mm fin thickness. This sequence holds for 
both pressure conditions. The rate of condensate flow 
around the tube was observed to be dependent on the fin 
thickness. The fastest flow pattern was observed for the 
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tube with 0.5 mm fin thickness, then for the tube with 0.75 
mm and finally for the tube with 1.0 mm fin thickness. This 
surprising observation can be explained if the fin tip 
effect is introduced; i.e., since the fin tip is rectangu- 
lar, there is a possibility that some of the condensate 
will fall into the interfin spacing and some will follow the 
circumferential path. If this condensate flow is realistic, 
it can be concluded that the thicker the fin, the more 
condensate will flow circumferentially along the fin tip. 

To further support this conclusion, numerous runs 
were made at a different fin spacing (1.0 mm) where the 
miniscus in the interfin area at the top part of the tube 
is visible (recall that the tubes with 0 . 5-mm fin spacing 
are completely flooded) . If the miniscus radius is different 
as the fin thickness changes, it can be concluded that the 
condensate flow pattern is influenced by the fin thickness. 
The tube with 0.5 mm fin thickness had a smaller miniscus 
radius than the tube with 1.0 and 1.5 mm fin thickness. 

This means that the previous conclusion may be correct; 
therefore, when the fin tip is wide and of course rectangu- 
lar, there will be circumferential drainage along the fin 
tip. On the other hand, for tubes with narrow fin tips, 
most of the condensate will fall into the interfin spacing, 
resulting in a very fast flow pattern and the creation of 
ripples. If this is the case, the thin fins should show 
an increased performance due to the convective action of 
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the condensate. However, the final performance is limited 
by the excessive condensate amount that is attached continu- 
ously on the tube. These may be the reasons why the tube 
with 0.5 mm fin thickness has better performance than the 
tube with 1.0 mm fin thickness when the fin spacing is 
0.5 mm. Between the tubes with 0.5 and 0.75 mm fin thick- 
nesses, there is no clear visible difference in the flow 
pattern; therefore, the last tube is superior because it 
presents less conduction resistance (it holds less conden- 
sate) , and it certainly has a higher fin efficiency. 

In general , it can be concluded that for fins of l._0 
mm height , the enhancement ratio was weakly dependent on fin 
thickness . The most significant reason is that the fin 
thickness has no influence on the retention problem as the 
fin spacing does . 

In order to check the validity of some of the 
previous conclusions, data were taken for one more tube. 

This tube had a fin spacing, thickness and height of 1.25 
mm, 1.25 mm and 1 mm, respectively, and it was made from the 
tube with 1.5 mm fin thickness whose fins were machined down 
by 0.25 mm. A tube with 1.25-mm fin spacing was expected 
to have a better heat-transfer performance than the tube of 
1.0 mm fin spacing since the condensing heat-transfer coeffi- 
cient increases up to a fin spacing of around 1.5 mm. 

However, the performance should be limited owing to the 
fact that the fin thickness is outside the optimum range 
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(0.75 to 1.0 mm) . Figures 5.28 and 5.29 show the relative 
plots of the condensing coefficient with heat flux for both 
pressure conditions, for three tubes whose fin geometry is 
shown in these figures. It is shown that the performance 
of the new tube is almost the same as with the tube with 
1.0 mm fin spacing and thickness. This means that what was 
gained by the spacing increment (from 1.0 to 1.25 mm) was 
lost by the fin thickness increase (from 1.0 to 1.25 mm) ; 
therefore the previous conclusions about the fin spacing and 
thickness effects seem to be reasonable. 

4 . Significance of Fin Height 

As shown in Table I, two fin heights (1.0 and 2.0 
mm) were examined during this thesis effort. The geometry 
of tubes tested is given in Table I. The heat-transfer 
performance of these tubes has already been shown with 
Figures 5.22 to 5.25, and for both sets of tubes, the 
optimum fin spacing was determined to be 1.5 mm. Figure 
5.30 shows the enhancement ratio for the set of tubes with 
2 . 0-mm fin height. It is shown that the maximum enhancement 
is about 5 and 6.5 for the vacuum and atmospheric runs, 
respectively. Figure 5.31 is a cross plot showing the 
enhancement ratio versus the fin spacing for both pressure 
conditions, with the fin height as a parameter. It is shown 
that the set of tubes with 2 mm fin height has better 
performance than the set of tubes with 1.0 mm fin height. 

The reasoning is two fold: 
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Figure 5.28. Variation of Condensing Coefficient with Heat Flux 

for Tubes with Fin Geometry Shown in the Figure 
(Vac. Runs) 



Atmospheric 
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Figure 5.29. Variation of Condensing Coefficient with Heat Flux 

for Tubes with Fin Geometry Shown in the Figure 
(Atm. Runs) 
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Figure 5.30. Cross Plot Showing Best Fin Spacing for Tubes with 

e = 2.0 mm and t = 1.0 mm, for both Pressure Conditions 
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Figure 5.31. Comparison of Enhancement Ratio Obtained on Sets 

of Tubes with e = 1.0 and e = 2 . 0 mm (Vac. Runs) 



1. Surface area increases with increasing fin height 
(see Figure 5.30) . 

2. The retention angle decreases and hence heat-transfer 
performance increases, with increasing fin height as pre- 
dicted by equation (1.12). This was also visually veri- 
fied during this thesis effort. 

Examining carefully the above mentioned Figure 5.31, 
it can be seen that the performance sequence is not the same 
for both sets of tubes. The tube with 1.0-mm fin spacing 
is the second in performance among the set of tubes with 
1.0-mm height. However, for the set of tubes with 2.0 mm 
height, the tube with 1.0-mm fin spacing is the third in 
performance while the tube of 2.0-mm fin spacing becomes 
second. This observation holds for both pressure conditions. 
This means that as the fin height increases, there is a 
tendency for the optimum spacing to mo\^ to a higher value. 
This tendency is mainly due to the increased amount of 
condensate which is retained between the fins at the top 
part of the tube, caused by the increased condensate 
generation. This fact was verified visually during the 
experiments, and Figure 5.32 shows a schematic of the con- 
densate appearance on those two tubes which differ only 
by the fin height. It is shown that the interfin space at 
the top part of the tube with 1.0-mm fin height 'holds' less 
condensate than the tube with 2.0-mm fin height. This obser- 
vation may reveal that the mechanism of thinning the film 
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Finned Tube with s = 1.5. t = 1.0 and 




113 



Figure 5.32. Condensate Profiles for Tubes with Fins 

of Different Height 






diminishes when the fin height increases. The pressure 
gradient in the condensate film between points A and B 
(see Figure 1.1) can be approximated as: 



dP \ 



AP 



AB 




e 



(5.13) 



Thus, the pressure gradient between points A and B decreases 
with increasing fin height (e) ; therefore, it should be 
expected that the tubes with 2.0-mm fin height would provide 
a less effective thinning effect. To prove this statement, 
one more cross plot was prepared by plotting the normalized 
enhancement ratio versus the fin spacing (i.e., separate 
the area addition effect). Figure 5.33 shows clearly that 
the set of tubes with 2.0-mm fin height is not improving the 
condensation enhancement mechanisms as effectively as 'the 
set of tubes with 1.0-mm fin height. 



increases, the portion of the fin that extends above the top 
flooded part of the tube increases (see Figure 5.32 and 
compare the extending fin part for the tubes with 1.0 and 
2.0 mm fin height); therefore, this part of the fin presents 
the highest heat-transfer performance since the condensate 
layer around the fin is extremely small. Of course, the 
performance of the tip is limited by the fin efficiency 
which decreases with increasing fin height; therefore, the 
set of tubes with 2.0-mm fin height has better performance 
mainly because of the excessive area addition. 



It is worth mentioning that as the fin height 
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Figure 5.33. Normalized Enhancement Ratio of Sets of Tubes 

with Fins of 1.0 and 2.0 mm Height 



Up to now, two ways to increase the external surface 
area of a tube were mentioned; i.e., either by reducing the 
fin thickness or by increasing the fin height. In the 
previous section it was shown that by changing the fin 
thickness there was not much enhancement, the reason being 
that the fin thickness does not improve the flooding problem 
and that the area addition happens in the regimes of the 
tube that are flooded. On the other hand, by_ increasing the 
fin height there is the benefit of area addition away from 
the flooded part of the tube resulting in a better heat- 
transfer performance . 

5 . Significance of Tube Diameter 

Yau et al. [Ref. 3] took data for filmwise conden- 
sation on copper finned tubes under near-atmospheric pres- 
sure. The fin shape was rectangular and the fins were 0.5 
mm in thickness, 1.0 mm in height, and the fin spacings were 
from 0.5 up to 20.5 mm (13 different fin spacings). The 
diameter of the tubes (at the root of fins) was 12.7 mm and 
the vapor velocity was 1.1 m/s. 

Figure 5.34 shows the enhancement ratio versus fin 
spacing for the data taken by Yau et al . and the data taken 
during this thesis effort. The only difference between the 
two sets of tubes was the root diameter, which was 12.7 and 
19 mm for Yau et al . and the present work, respectively. 

It can be seen that there is a similarity in the enhancement 
ratio curves and for both sets of tubes the best fin spacing 



116 




ID 



^r 



co 



CM 



O 



E 

E 

\ 

CD 

C 

•H 

0 

0 

CL 

CO 

c 

T"J 

LL 



(0 LD -t 



CO CM 



o 



OT^Dy ^UBUJ0OUDL)U3 



117 



Figure 5.34. Cross Plots Showing Best Fin Spacing for 



found was 1.5 mm. However, there is a significant differ- 
ence in the enhancement ratio, which can be explained by 
the following reasons: 

1. The retention angle increases as the tube diameter 
decreases, as it is predicted by equation (1.12), resulting 
in a lower condensing coefficient. 

2. The smooth-tube performance is predicted by the 
Nusselt equation, and it gives that a tube with smaller 
diameter will have higher condensing coefficient (recall 
that the tube diameter is in the denominator) . 

Therefore, when the enhancement ratio is computed, 
the tube with smaller diameter will have the numerator small 
and the denominator large, resulting in a smaller ratio 
than that which was obtained during this thesis effort. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

1. For the best fin spacing of 1.5 mm, a fin thickness 
of 1.0 mm and a fin height of 2.0 mm, heat- transfer enhance- 
ments around 4 to 5 were realized for the vacuum runs and 
around 6 for the atmospheric runs, while the area increment 
was only about 3 times the smooth- tube area. The additional 
enhancement is due to the thinning effect; therefore, any 
theoretical model that neglects this effect may have severe 
deviations from the experimental data. 

2. The 1.5 mm fin spacing appears to be the optimum 
among the fin spacings tested (0.5, 1.0, 1.5, 2.0 and 4.0 mm) 
for both vacuum and near-atmospheric pressure conditions. 

3. The condensation enhancement does not strongly depend 
on the fin thickness. However, the tubes with fin thickness 
around 0.75 to 1.0 mm performed better than the tubes with 
0.5 and 1.5 mm fin thicknesses. 

4. As the fin height increased from 1.0 to 2.0 mm, 

the condensing coefficient increased only by about 20 percent, 
despite an area increase of about 50 percent. This lower- 
than-expected performance is believed to be owing to the 
poorer thinning effect of the condensate film with increas- 
ing fin height. 

5. As the tube diameter increases, the retention angle 
decreases resulting in a better heat-transf er performance; 
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therefore, the tubes with a root diameter of 19.0 mm showed 
a superior enhancement over the tubes with a 12.7 mm root 
diameter . 

6. The consistency of equation (1.12) that predicts the 
retention angle was visually verified throughout all the 
experiments. Certainly, it must be a very useful equation 
to predict the condensate retention angle when modelling 
steam condensation on finned tubes. 

7. The leading coefficient of the Sieder-Tate-type 
equation used- to determine the water-side resistance with 
the instrumented tube was 0.0635, while the modified Wilson 
plot technique gave the same coefficient as 0.071. 

8. The test apparatus successfully operates under any 
pressure condition. An extremely tight apparatus is 
achieved resulting in minimization of the non-condensing 
gas effect. 

9. A tube-cleaning and preparation procedure has been 
established to ensure filmwise condensation on copper tubes. 

B. RECOMMENDATIONS 

1. Calibrate the existing pressure transducer in order 
to replace the U-tube, mercury-in-glass manometer. This 
will make each run faster and more automatic (less data to 
be entered manually) . 

2. Install the relief valve on the test condenser. 

3. Take data with different shaped fins (at the tip as 
well as at the root of the fin) to investigate the effect 
of fin geometry. 
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4. Take data with tubes that are fabricated with 



different thermal conductivity materials to investigate the 
effect of the material conductivity on the condensing 
heat-transfer coefficient. 

5. Attach drainage strips on the existing tubes to 
investigate the effect of the drainage strips on steam 
condensation enhancement. 

6. Take data to study the effect of vapor shear on 
finned tube enhancement. For the vacuum runs, it is recom- 
mended that vapor velocities of 1.0, 2.0, 3.0 and 4.0 m/s 
be used and for the atmospheric runs the vapor velocities 
of 0.7 and 1.0 m/s. Higher vapor velocities would be 
desirable if the apparatus could be modified. 

7. Modify the existing apparatus to study the conden- 
sate inundation effect and take data to investigate if the 
enhancement and the optimum fin spacing are affected by 
inundation . 

8. Take data on finned tubes that are coated with a 
dropwise promoter. Since the presence of the fins will 
result in a more effective bonding of the promoter, this 
combination may result in superior enhancement and promoter 
endurance. Also, the condensate may not be retained on 
dropwise-coated finned tubes with sufficient fin spacing. 

9. Modify equation (1.14) of Owen et al. [Ref. 8] 

to include the surface-tension effect on the unflooded part 
of the tube. By including this effect, the Owen prediction 
for the steam data will become less conservative. 
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APPENDIX A 



PREPARATION OF CONDENSER TUBES 



Before every run, the tube was prepared by following 
the steps below. The operator must wear gloves to avoid 
contamination of the tube as well as to avoid skin irritation 
that may be caused by the cleaning solution. 

1. Prepare a solution of equal parts of sodium hydroxide 
and methyl alcohol. Heat the solution until it becomes 
watery. Using a toothbrush, apply a coating of this mixture 
on the tube. 

2. Place the coated tube in a steam bath for about one 
hour . 

3. Rinse the tube with tap water to avoid contamination 
of the apparatus with the. solution. 

4. Clean the external surface of the tube with a brush 
using a weak solution of tap water and a detergent (Spark- 
leen) . This will ensure that all the sodium hydroxide 
deposits are removed from the corners of the fins. 

5. Using a brush, clean the inside surface of the tube. 

It is essential to remove any fouling. 

6. Install the tube as soon as possible to avoid any 
system contamination from the environment. 
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APPENDIX B 



SYSTEM START-UP AND SHUT-DOWN PROCEDURES 

To start the system, the following procedure is 
recommended : 

1. Ensure that the boiler water level in the boiler is 
at least six inches above the upper level of the heaters. 

2. Close the venting valve and open the valve on the 
purge line (see Figure 2.1). 

3. Ensure that the vacuum pump is filled with oil at the 
right level. 

4. Turn on the vacuum pump. 

5. If the run is for low pressure then turn on the boiler 
power when the system pressure approximates the desired 
pressure. Avoid boiling if the system pressure is greater 
by 50 mmHg than the desired pressure. Boiling at an ele- 
vated pressure causes the apparatus walls to reach a 
temperature higher than the desired saturation temperature. 
If this happens, a considerable time interval (as much as 30 
minutes for a temperature difference of 10 K) will be needed 
to avoid superheated steam. Note that when the system walls 
are at an elevated temperature, the steam gains heat, thus 
causing superheating. The power input is controlled by a 
potentiometer until the voltmeter reads 90 units (note that 
each unit represents 2.2 Volts). This setting yields a 2.0 
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m/s vapor velocity at a system pressure of 85 mmHg . If the 
run is an atmospheric one, then, turn on the power right 
after the vacuum pump is energized. For a vapor velocity 
of 1.0 m/s, the voltmeter should be set at its maximum 
which is 170 units or (170.0 x 2,2) Volts. 

6. Open the valve to fill the tap water sump. 

7. Circulate tap water through the auxiliary condenser. 

8. Check that the water sump is filled and then turn on 
the pumps that circulate water through the test bue. 

9. When the steam reaches the tube (as can be seen from 
the viewport) , adjust the water flowrate through the auxili- 
ary condenser in order to achieve the desired pressure. 

10. Energize the data acquisition system, the computer, 
the printer and insert in the driver of the computer the 
disk which contains the code "DRP4." Finally, load "DRP4" 
and follow the interactive program to gather and analyze 
data . 

To secure the system: 

1. Turn off the input power to the boiler. 

2. Secure the water supply to the auxiliary condenser. 

3. Open the venting valve to relieve the pressure in the 
apparatus . 

4. Turn off the printer, the data acquisition system 
and the computer. 

5. Secure the pumps that circulate tap water through 
the test tube. 
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APPENDIX C 



RAW DATA OF UNINSTRUMENTED AND INSTRUMENTED TUBE 



The following two tables contain the raw data taken for 
the uninstrumented and instrumented tubes under vacuum and 
near-atmospheric conditions. 



TABLE V 

Raw Data of Uninstrumented Smooth Tube 



File Name 

Pressure Condition 
Vapor Velocity 



Vw 


Tin 


m/s 


°C 


4.47 


18.58 


3.93 


18.63 


3.38 


18.70 


3.10 


18.71 


2.71 


18.81 


1.99 


18.99 


1.48 


19 . 19 


1.15 


19.46 



File Name 

Pressure Condition 
Vapor Velocity 



Vw 


Tin 


m/s 


°C 


4.49 


18.84 


3.93 


18.91 


3.38 


18.99 


3.10 


19.00 


2.71 


19.11 


1.99 


19.32 


1.48 


19 . 56 


1.15 


19.91 



WSIV14 

Vacuum (85 mmHg) 
1.0 m/s 



Tout 


Ts 


°C 


°C 


19.37 


49.45 


19.48 


49 . 59 


19 .65 


49.46 


19 .72 


49.76 


19.92 


49.41 


20 . 34 


49 . 84 


20 . 87 


49.49 


21.46 


49 . 53 


WSIA16 




Atmospheric 




0 . 8 m/s 




Tout 


Ts 


°C 


°C 


20 . 65 


99 . 81 


20.94 


99.96 


21.31 


100.11 


21.48 


99.93 


21.85 


100.31 


22 . 84 


100 . 26 


23.92 


99.98 


25.20 


99.97 
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Data of Instrumented 
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APPENDIX D 



RAW DATA OF TUBES TESTED 



The following pages contain tables with the raw data 
taken for all tubes tested under vacuum and near-atmospheric 
conditions . 
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APPENDIX E 



ERROR ANALYSIS 



The numerical data taken during an experiment are used 
together with theoretical formulations in order to obtain 
results about the objective of the experiment. Any measure- 
ment has an uncertainty which is a function of the device 
accuracy as well as the experience of the operator. 
Therefore, the final result may be severely distorted due 
to the error propagation during the computational steps. 

In cases where the final result has a large uncertainty, 
it may be unwise to draw any conclusions. 

The probable uncertainty in a computation can be 
determined from the following equation, as taken from 
[Ref. 19] : 



w 



R 



= [( 



3R 

9x. 



Wj^) 



+ ( 



9R 

9x, 



w 2 ) 



+ 



+ ( 



9R_ 

9x„ 



n 



w ) ] 

n 



1/2 



(E.l) 



where : 



R is the result 

x^,X 2 /...,x are the measured (independent) 
variables 

w l, w 2,.**, w n are t * le uncertainties in the 
measured (independent) variables. 

The above equation, together with the uncertainty of 
each measurement, was used to determine the possible error 
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of the outside heat- transfer coefficient. As discussed in 



Chapter V, the outside resistance was found by subtracting 
the water-side and tube-wall resistances from the overall 
resistance. This means that the error in the condensing 
coefficient is mainly a function of the overall and the 
water-side resistances. In such cases, if the water-side 
resistance has a significant error it will be unavoidably 
reflected in the condensing heat-transfer coefficient. 
Further, if the value of the water-side resistance is much 
larger than the condensing resistance, then the error that 
will appear in the condensing coefficient will be highly 
amplified. The steam-side resistance is generally smaller, 
especially on finned tubes, then the water-side resistance 
(recall that the condensation process is a phase change 
phenomenon where the heat-transfer process is very effec- 
tive) . Therefore, to 'protect' the final result of the 
condensing coefficient from an excessive error, the water- 
side and condensing resistances should be as nearly balanced 
as possible. The water-side resistance can be reduced 
either by increasing the coolant velocity or by internal 
enhancement. Both ways promote turbulence and reduce the 
boundary layer thickness on the inside wall surface (recall 
that the coiled insert was used exactly for this purpose) . 
However, since the water-side temperature rise decreases 
with increase in water velocity, the water velocity must not 
be too excessive. Since the water-side temperature rise is 
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the most crucial measurement, a minimum value must be 
assigned to minimize experimental errors. For this reason, 
the internal enhancement is preferable than the increase 
in water velocity to obtain a lower thermal resistance on 
the waterside. 

Based on the above discussion, it should be expected 
that as the cooling water flow rate decreases (i.e., less 
turbulence) , the error in the outside heat-transf er coeffi- 
cient should increase. Figures E.l and E.2 show the error 
bands for three cooling water flow rates for the tube that 
showed the best heat- transfer performance (1.0 mm fin thick- 
ness, 1.5 mm fin spacing and 2.0 mm fin height) for both 
pressure conditions. It is shown that indeed the error 
gets larger as the cooling water flow rate (i.e., the heat 
flux) decreases. Also notice that the error bands shown 
with the atmospheric run are smaller than the corresponding 
value for the vacuum runs. This improvement can be ex- 
plained by the increased water-side temperature rise. 

Note that the relative error decreases as the water-side 
temperature rise increases. 

Based on the foregoing discussion, it is clear why the 

reference heat fluxes for computing the enhancement ratio 

of the cross plots in Chapter V were selected as high as 

2 

possible (0.35 and 1.0 MW/m for vacuum and atmospheric 
runs, respectively). 

The uncertainty of the major variables that govern the 
condensation process (cooling water velocity, Reynolds 
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Figure E.l. Error Bands for a Tube Tested Under Vacuum Conditions 
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Figure E.2. Error Bands for a Tube Tested Under Atmospheric Conditions 



number, heat flux, LMTD, water-side and outside heat- 
transfer coefficients) are given as percentages in the 
following computer outputs of the program "ERROR. " At 
the end of this appendix, the listing of the program 
"ERROR" is provided for the Error Analysis. 
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OA T A FOR the ERROR ANALYSIS: 



File Name: F25E2VM6 

Pressure Condition: VACUUM !i 

Coding Water Temperature In 
Cooling Water Temperature Out 
Steam Temperature 
Water Floy Rate (7) 

Water Veioci ty 



k Pa 
22.72 
24.07 
48.35 
80.00 
4.39 



( D e g C ) 
(Deg C ) 
(Deg C ) 



( m/s ) 



L'-'ROR ANALYSIS: 



7 


Error 


i n 


Mass 


Floy Rate 


= 


1 . 03 


7 


Error 


l n 


Reyno 


Ids Number 


= 


i .24 


/• 


Error 


l n 


Heat 


F 1 UX 


= 


0 . 79 


X 


E r r o r 


i n 


LMTD 




= 


2.22 


7. 


Er ro r 


in 


iJo 




= 


2.38 


7. 


Error 


i n 


Hi 




= 


1 .32 


V 

/> 


Error 


in 


Ho 




= 


8 . 36 



DATA FOR THE ERROR ANALYSIS: 



tile Name: F25E2VI46 

Pressure Condition: VACUUM I I kPa 

r « 1 I i — j- _ t j t _ - i o 



Cooling Water Temperature 


In 


23.02 


( Deg 


C> 


Coding Water Temperature 


Out = 


24 . 79 


( Deg 


C) 


Steam Temperature 


= 


48.20 


(Deg 


c> 


Water Floy Rate (7.) 


= 


45.00 






Water Velocity 


= 


2.51 


( m/s 


) 



ERROR ANALYSIS: 






7. Err o r 


i n 


Mass Floy Rate 


= 


i .80 


7 Error 


in 


Reynolds Number 




Q p 


7 Error 


i n 


Heat Flux 


= 


0.78 


7 Error 


m 


LMTD 


= 


1 .63 


7 Err o r 


l n 


U o 


- 


1.3! 


7 Error 


in 


Hi 


= 


1 . 78 


7. Error 


i n 


Ho 


= 


* 0.58 
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DATA FOR THE ERROR ANALYSIS: 



Name : 

•sure Condition 



F25E2V146 
VACUUM 1 I 



r i ie 
Pr e 

Cooling Hater Temperature In 
Cooling Hater Temperature Out = 
Steam Temperature 
Hater Flow Rate <X) 

Water Velocity 



Ie Pa 



23 , 


.53 


< Deg 


C> 


26 , 


, ! ! 


( Deg 


C) 


48 , 


.53 


( Deg 


C) 


20, 


. 0 0 






1 , 


.16 


< m/s 


) 



ERROR ANALYSIS: 



'•/ 
* > 


Error 


i n 


Mass Flow Rate 


= 


O O Q 
. J • <. > I* 


X 


Error 


in 


Reynolds Number 


= 


3.96 


X 


Error 


in 


Heat Flux 


= 


1 . 30 


x 


Error 


in 


LMTD 


= 


1.10 


X 


Error 


i n 


Uo 


- 


* .70 


y 


Error 


i n 


Hi 


= 


O O Q 


X 


L r r o r 


j.n 


Ho 


- 


36.25 



DATA FOR THE ERROR ANALYSIS: 

File Name: F25E2A128 

Frees ure Condition: ATMOSPHERIC 101 ’<Pa 



Cooling Hater Temperature 


In - 


22.90 


( Deg 


C ) 


Cooling Hater Temperature 


Out = 


26.96 


< Deg 


C) 


Steam 


Tempe nature 


= 


99.37 


( Deg 


C> 


Hater 


Flow Rate (X) 


= 


80 . 00 






Water 


Velocity 


= 


4.39 


( m/s 


> 



ERROR ANALYSIS: 






X Error 


i n 


Mass Flow Rate 


_ 


1 .03 


X Error 


i n 


Reynolds Number 


= 


1 .25 


X Error 


i n 


Heat Flux 




0.42 


X Error 


i n 


LMTD 


= 


0.71 


X Error 


i n 


Uo 


= 


0.83 


X Error 


in 


Hi 


= 


1 .33 


X E r r o r 


in 


Ho 


- 


3 . 5 I 
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DATA FOR THE ERROR ANALYSIS 



File Name : F25E2A I 28 

Pressure Condition: ATMOSPHERIC 

Cooling Hater Temperature In 
Cooling Hater Temperature Out = 
Steam Temperature 
Water Flow Rate <7») 

Water Velocity 



1 0 1 k Pa 
23.27 (Deg C> 
28.87 (Deg C) 
98.59 (Deg C> 
45.00 

2.51 (m/s) 



ERROR ANALYSIS: 



7 

/ » 


Error 


in 


Mass Flow Rate 


= 


> .30 


7 


Error 


i n 


Reynolds Number 


= 


1 ,95 


V 

/> 


Error 


i n 


Heat Flux 


= 


0 . 6 ! 


7. 


Error 


in 


LMTD 


= 


0 . 5 1 


% 


Error 


i n 


Uo 


= 


0.80 


7 


Error 


l n 


Hi 


- 


1 .79 


% 


Error 


i n 


Ho 


= 


6 . 34 



DATA FOR THE ERROR ANALYSIS: 

File Name: F25E2AI28 

Pressure Condition: ATMOSPHERIC 

Cooling Water Temperature In 
Cooling Water T emperature Out = 
Steam Temperature 
Water Flow Rate (7.) 

Water Velocity 



23.96 (Deg C> 
32.30 (Deg C) 
99.43 (Deg C> 
20.00 

1.16 ( m/s > 



ERROR ANALYSIS: 



7. E r r o r 


in 


Mass Flow Rate 


= 


3.90 


7 Error 


in 


Reynolds Number 


= 


3.97 


7. Error 


i n 


Heat F lux 


- 


1 .25 


7 Error 


in 


LMTD 


= 


0.34 


7. Error 


m 


Uo 


= 


1 .30 


7 Error 


in 


H i 


= 


3.30 


7 Error 


i n 


Ho 


= 


27.01 
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! 0 0 0 ! 

1010 ! 
1020 ! 
1030 
1040 
1050 
1080 
1070 
1080 
1090 
1100 
1 1 to 
1120 
1130 

1 14(5 
1150 
1 180 
1170 
1180 
Tci 
1 190 
1200 
1210 
1220 
Tco 
1230 
1240 
1250 
Ts 
1260 
1270 
1280 
1290 
1300! 
1310 
1320 
1330 
1340 
1350 
1360! 
1370 
1380 
1390 
1400 
1410 
1420! 
1430 
1440' 
w 

1450! 

1460 

1470 

1480! 

1490 

1500! 

1510 

1520 

1530 

1540 



FILE NAME : ERROR 
REVISED : JUNE 27. 1984 



PRINT USING ”i OX. '••’DATA FOR THE ERROR ANALYSIS: 

PRINT 

BEEP 

INPUT "ENTER FILE NAME". FiieS 

PRINT USING "1 OX . ""File Name: "" . 12A" ;Fi le$ 

BEEP 

INPUT "ENTER PRESSURE CONDITION < 1 =VACUUM .2=ATM0SPHERIC> ".Prc 

IF P rc = 1 THEN 

BEEP 

PRINT USING "10X. ""Pressure Condition: VACUUM 11 kPa 

ELSE 

PRINT USING "10X. ""Pressure Condition: ATMOSPHERIC 101 kPa 

END IF 
BEEP 

INPUT "ENTER COOLING WATER TEMPERATURE IN". Tci 

PRINT USING "1 OX .""Cooling Water Temperature In = "".3D.2D,"" (Deg C) 
BEEP 

INPUT "ENTER COOLING WATER TEMPERATURE DIFFERENCE" , Dtu 
Tco=Tci+Dtu 

PRINT USING " 1 0X . ""Coo 1 ing Water Temperature Out = "".3D. 2D,"” (Deg C> 



BEEP 

INPUT "ENTER STEAM TEMPERATURE" . Ts 

PRINT USING "10X, ""Steam Temperature = "".3D.2D,"" (Deg C) : 



BEEP 

INPUT "ENTER COOLING WATER FLOW RATE 7",Fm 

PRINT USING "1 OX. ""Water Flow Rate (7.) = "”,3D.2D";Fm 

T=(Tci+Tco)/2 ! FILM TEMPERATURE 

UNCERTAINTY IN THE COOLING WATER 

Df m= . 1 0 ! ERROR IN FLOW METER READING 

Drho= 1 ! ERROR IN WATER DENSITY 

Dmf=5 . 00049E-3+6 .9861 937E-3*Dfm ! ERROR IN MASS FLOW RATE 
Rho=FNRho(T) ! WATER DENSITY 

Mf=1 . 04805E-2+6 . 80932E- 3»Fm ! MASS FLOW RATE OF COOLING WATER 

CORRECT MF FOR THE TEMPERATURE EFFECT 

Mf =Mf*( 1 .0365-1 .96644E-3*Tci+5.252E-6*Tci'2>/. 995434 

D i = .0127 ! TUBE INSIDE DIAMETER 

A i =( P I»D i ~2 ) /4 ! TUBE INSIDE CROSS SECTION AREA 

Ddi=. 000025 

Dai=<PI*Ddi 2>/4 ! ERROR OF INSIDE TUBE CROSS AREA 
COMPUTE THE WATER VELOCITY 
Vu=Mf / (Rho*Ai ) ! WATER VELOCITY 

PRINT USING "1 OX , ""Water Velocity = , Z . DD ."" (m/s) 



V 



CORRECT OUTLET WATER TEMP. FOR THE MIXING CHAMBER EFFECT 
Tco=Tco-.004*Vw'2 ! COOLING- WATER OUTLET TEMP. 

T=( Tc i + Tco) » . 5 ! FILM TEMPERATURE 
COMPUTE THE ERROR IN WATER VELOCITY 
Dvw=Vu*SQR ((Dmf/Mf) '2 + ( Drho/Rho > '2 + ( Dai/A i ) '2 ) 

UNCERTAINTY IN THE REYNOLDS NUMBER 
Mu=FNMu< T) ! WATER VISCOSITY 
Dmu=6.E-6 ! ERROR OF WATER VISCOSITY 
Re= ( Rho»Vw*D i ) /Mu 

Dre=Re*SQR ( (Drho/Rho) '2 + ( Dvu/Vu ) 2 + <Ddi/Di) '2 + ( Dmu/Mu ) 2 ) 
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1550! UNCERTAINTY IN THE HEAT TRANSFERED 
1560 Cpw=FNCpu ( T ) 

1570 Q=Mf * < Tco- Tci > *Cpu 
1580 Dtco=.02 
1590 Dtci=.02 
1600 Dcpu=8 

1610 Dq=Q*SQR ( ( Dmf /Mf ) ' 2+< ( D tco/ (Tco-Tci))) 2 + ((Dtci/(Tco-Tci)))'2 + <Dcpii»/Cpw) '2 
) 

1620! UNCERTAINTY IN THE HEAT FLUX 
1630 Dl= . 0005 ! ERROR IN TUBE LENGHT 
1640 Do-. 01905 ! SMOOTH TUBE OUTER DIAMETER 
1650 Ddo=. 000025 

1660 L=. 13335 ! FINNED TUBE LENGTH 

1670 Q P =Q/ <PI*Do*L) ! HEAT FLUX 

1680 Dqp=( Qp/PI ) *$QR ( ( Dq/Q) ~2 + < Ddo/Do )''2 + (DI/L)'‘2) 

1690 Dts = . 02 ! ERROR IN TEMPERATURE MEASUREMENT 

17(10 Lmtd=(Tco-Tci) /L0G( <Ts-Tci)/(Ts-Tco)) 

1710 Uo=Qp/Lmtd ! OVERALL HEAT TRANSFER COEF. 

1720 A1=Dts*(Tci-Tco)/((Ts-TcL)*<Ts-Tco) *L0G< (Ts-Tci)/(Ts-Tco)>) 

1 730 A2=Dtci/((Ts-Tci)*L0G((T3-Tcx)/<Ts-Tco) )) 

1740 A3=Dtco/( (Ts-Tco)*LOG<<Ts-Tci)/(Ts-Tco) ) ) 

1 750 DImtd=Lmtd*SQR< A1 ~2+A2~2+A3'2> 

1 760 Duo=Uo*S0R< (Dqp/Qp) "2+< Dlmtd/Lmtd) '2) 

1770 M=Mu 

1780 T 1 =< T+273 . 15)/273. 15 
1790 Kw=FNKw ( T 1 ) 

1800 Ci = . 0635 ! SIEDER-TATE COEFFICIENT 

1810 Ac=26 . 4 ! INTERSCEPT FROM SIEDER PROGRAM 

1820 Kc=385 ! WALL THERMAL CONDUCTIVITY 

1830 LI =.060325 ! LENGTH OF UNFINNED LEFT PART OF TUBE 

1840 L2=. 034925 ! LENGTH OF UNFINNED RIGHT PART OF TUBE 

1850 Dr-. 015875 

1360 Pr=Cpu*Mu/Ku) 

1870 Muu) = FNMuu( T ) 

1380! UNCERTAINTY OF INSIDE HEAT-TRANSFER COEFF . 

1890 Cf-1. 

1900 H i = ( Kw/D i > -*< C i*Re . 8*Pr ' .333*Cf +Ac) 

1910 Dti=Q/(PI»Di*(L+L1»Fe1 +L2*Fe2>*Hi) 

1920 Cf c = < Muoj/FNMuu ( T+D 1 1 ) ) ' . 1 4 

1930 IF ABS( (Cfc-Cf )/Cfc)> .01 THEN 

1940 Cf-(Cf+Cfc)*.5 

1950 GOTO 1900 

1960 END IF 

1970 PI =PI*<Di+Do) 

1980 Bl-<Do-Di)*PI*(Di+Do)*.5 
1990 M1=(Hi*P1 /<Kc*91 > >' .5 

2000 P2=PI*< Di + Dr > 

2010 B2=< Dr-Di )*PI*(Di+Dr)*.5 

2020 M2=<Hi»P2/<Kc~B2> ) ' .5 

2030 Fel =FNTanh(M1 *L1 >/<M1 *L1 ) 

2040 Fe2=FNTanh(M2*L2)/(M2*L2) 

2050 Dtc=Q/(PI*Di*(L+L1*Fe1 +L2-Fe2)-Hi> 

2060 IF ABS<<Dtc-Dti)/Dtc)>.01 THEN 1900 

2070 Die u= . 001 0 ! ERROR IN WATER THERMAL CONDUCTIVITY 

2080 Dci= .0005 ! ERROR IN SIEDER-TATE COEFFICIENT 

2090 Dpr= . 05 ! ERROR IN PRANDTL NUMBER 

2100 Dcf =8 .E-6 

2110 A4= . 1 4*Dcf /Cf 

2120 Dh i=H i *SQR < ( Dku/Kw) '2 + < Dd i/D i ) 2 + ( . S^Dre/Re ) '2 + < .333*Dpr/Pr) 2 + <Dci/Ci) '2 + 
A4 > 

2130! UNCERTAINTY OF OUTSIDE HEAT-TRANSFER COEFF. 
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2140 Dkc = 1 ! ERROR IN TUBE THERMAL CONDUCTIVITY 

2150 Ru=Do*LGG< Do/D 1 ) /< 2*Kc) ! WALL RESISTANCE 

2160 Ho-1 / < <1/Uo)-<Do*L/<Di*(L+L1*Fe1 +L2*Fe2)*Hi ) ) -Ru> 

2170 Dru=Ru*S'QR ( ( Ddo/Do ) "2+< Dkc/Kc) '2 + < Ddo/ ( Do*L0G < Do/D i ) ) > 2 + < Dd i / ( D i*L0G < Do/D 
i) ) >"2> 



2t80 A5= I /Uo-Rui- < Do/ < D i*Hi ) > 

2190 A6=Duo/<Uo'2*A5> 

2200 A7 =Dru/A5 

2210 A8=< <Do/(Di*Hi> )*<Dhi/Hi) )/A5 

2220 PRINT 

2230 Dho=Ho»SQR(A6~2+A7~2+A8~2) 

2240! CALCULATE THE 7. ERROR IN Ho 
2250 Prho=Dho*1 00 /Ho 

2260! CALCULATE THE V. ERROR IN REYNOLDS NUMBER 
22 7j) Prre=Dre*1 00/Re 

2280! CALCULATE THE X ERROR IN MASS FLOW RATE 
2290 Prmf =Dmf *1 00/Mf 

2300! CALCULATE THE X ERROR IN HEAT TRANSFER 
2310 Prqp=Dqp*1 00/Qp 
2320! CALCULATE THE '/. ERROR IN LMTD 
2330 Pr lmtd=Dlmtd*1 00/Lmtd 

2340! CALCULATE THE X ERROR IN OVERALL HEAT TRANSFER C0EF. 
2350 Pruo=Duo*1 00/IJo 

2360! CALCULATE THE X ERROR IN INSIDE HEAT TRANSFER C0EFF. 



2370 Prhi=Dhi*1 00/Hi 
2380 PRINT 

2390 PRINT USING "1 OX', “"ERROR ANALYSIS: 

2400 PRINT 

2410 PRINT USING "10X. M "X Error in Mass Flow Rate = "" .2D . 2D . " : Prmf 

2420 PRINT USING "10X.""X Error in Reynolds Number = "" .2D . 2D . " :Prre 

2430 PRINT USING ,, 10X. ,,M X Error in Heat Flux = . Z . 2D . " ;Prqp 

2440 PRINT USING "10X.""X Error in LMTD = "" , Z . 2D" :Pr Imtd 

2450 PRINT USING "10X, M,, 7. Error in Uo = ,M * . Z . 2D . " :Pruo 

2460 PRINT USING "10X. M,, 7. Error in Hi = "" ,2D .2D . " :Prhi 

2470 PRINT USING ,, 10X, ,,, 7. Error in Ho = "" . 2D . 2D . " : Prho 



2480 END 

2490 DEF FNMuu(T) 

2500 A=247 . 8/ < T + 1 33 . 1 5) 

2510 Muu=2.4E-5*l0'A 

2520 RETURN Muu. 

2530 FNEND 

2540 DEF FNTanh < X ) 

2550 P=EXP(X) 

2560 0=EXP< -X) 

2570 Tanh= ( P-Q ) / (P+Q) 

2580 RETURN Tanh 
2590 FNEND 
2600 DEF FNKu< T 1 ) 

2610 Kuj=- .9224 7 +T1 *( 2 . 8395 -T 1 ♦< I .8007- T1 . 52577-. 07344*T1 )) ) 

2620 RETURN Ku 

2630 FNEND 

2640 DEF FNMu(T) 

2650 A=247 . 8/ < T+ 1 33 . 15) 

2660 Muj= 2 . 4E-5*1 0' A 
2670 RETURN Mu 
2680 FNEND 
2690 DEF FNRho< T) 

2700 Rho=999. 52946+ T*< . 0 1 269- T*(5 . 48251 3E-3-T* I .234147E-5) > 

2710 RETURN Rho 
2720 FNEND 
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J OJ c 



2730 



750 



^760 



DEF FNCpw ( T ) 

Cd>j=<4.21 1 20358-T*< 2 . 26826E-3-T*< 4 . 4236 1 E- 5+2 . 7 1 428E-7*T) » ) * 1 000 

RETURN Cpu 

FNEND 
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